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Abstract 

A unified theory of orthogonal polynomials of a discrete variable is presented 
through the eigenvalue problem of hermitian matrices of finite or infinite dimensions. 
It can be considered as a matrix version of exactly solvable Schrodinger equations. 
The hermitian matrices (factorisable Hamiltonians) are real symmetric tri-diagonal 
(Jacobi) matrices corresponding to second order difference equations. By solving the 
eigenvalue problem in two different ways, the duality relation of the eigenpolynomials 
and their dual polynomials is explicitly established. Through the techniques of exact 
Heisenberg operator solution and shape invariance, various quantities, the two types 
of eigenvalues (the eigenvalues and the sinusoidal coordinates), the coefficients of the 
three term recurrence, the normalisation measures and the normalisation constants etc. 
are determined explicitly. 



1 Introduction 

Due to the long history and rich and diverse applications, there are many different ways to 
introduce orthogonal polynomials [1]. Among many types and kinds of orthogonal polynomi- 
als, we will focus in this paper on the so-called orthogonal polynomials of a discrete variable 
[21, [3l m [5] , whose orthogonality measures are concentrated on discrete points, either finite 
or infinite in number. We will present a unified theory of these polynomials based on the 



eigenvalue problem (12. ip of a special class of hermitian matrices, which are real symmetric 
tri-diagonal (Jacobi) matrices (12 ■4p . (12. lip . Since the spectrum of a Jacobi matrix is simple, 
the orthogonality of eigenvectors is guaranteed. The eigenvalue problem of a Jacobi matrix 
can be considered as a difference equation version of the Schrodinger equation for one degree 
of freedom system, the most basic equation of quantum mechanics Various classical or- 
thogonal polynomials (the Hermite, Laguerre, Jacobi and their restrictions) have appeared 
as the eigenfunctions of exactly solvable quantum mechanics [6]. 

Roughly speaking, our line of arguments is a deformation (discretisation) of the main 
trend of the twentieth century mathematical physics; to pursue the parallelism between the 
matrix eigenvalue problem and the ordinary differential equations of the Sturm-Liouville 
type, with differential equations replaced by difference equations. 

We apply many ideas and methods for solving Schrodinger equations: Crum's theorem 
|7j, the factorisation method [8] (or the so-called supersymmetric quantum mechanics [9J), 
the method of exact Heisenberg operator solutions and the creation/annihilation operators 
[Tot [TT] together with symmetries (shape invariance (14.21) [T2| [13] and closure relation (14.310 
[To] ) to elucidate the universal structure of the orthogonal polynomials of a discrete variable. 
The Jacobi matrix eigenvalue problem Tiip = Sip (12. ip can be solved in a different way, 
through the three term recurrence relations for the dual polynomial in S. Combining the 
duality, which can be stated at many different levels, with the above mentioned solution 
techniques, various quantities of the orthogonal polynomials, the coefficients of the three 
term recurrence, the orthogonality measure and the normalisation constants, the difference 
equations for the dual polynomials, etc. are given explicitly in an elementary manner. We 
stress that duality is established equally for the finite and infinite dimensional cases. All the 
examples discussed in this paper are taken from the review of Koekoek and Swarttouw [2] 
and they are known to satisfy the closure relation (14.311) . The orthogonal polynomials with 
the Jackson integral type measure, e.g. the Big g- Jacobi polynomial, etc. will need different 
treatment and will be discussed elsewhere. As an inverse step to characterise the orthogonal 
polynomials of a discrete variable, we solve the closure relation algebraically to determine 
the possible forms of the Jacobi matrices. 

It should be stressed that the present duality has far richer contents than the well-known 
duality of the column and row eigenvectors. That is, in the (hermitian) eigenvalues problems, 
the orthogonality of the complete set of column eigenvectors automatically implies that of the 
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row eigenvectors. In the present case, the column eigenvectors are polynomials in sinusoidal 
coordinate fl2.28p r]{x) and the row eigenvectors are polynomials in the eigenvalues (12.281) 
S{n), and the duality is formulated as the equality of these two polynomials on the integer 
lattice points fl3.10p . 

The formulation of the unified theory of orthogonal polynomials based on the eigenvalue 
problems of Jacobi matrices is new. All the quantities and formulas are derived from the two 
input functions B{x) and D{x) (12. 4p . (12.51) . However, a good part of the explicit formulas 
for specific polynomials explored in section five is already known individually in the existing 
theories of orthogonal polynomials based on various settings. 

This paper is organised as follows. In section two, the general form of the hermitian 
matrices (or the Hamiltonians) is given and the solution procedure to arrive at the orthogonal 
polynomials is outlined. They are polynomials in the sinusoidal coordinate (I2.28P which plays 
an important role in exactly solvable quantum mechanics [TOl[ll]. The orthogonality measure 
is given by the solution (f)o{x) of the factorised equation (12.170 . The (po is the ground state 
wavefunction. In section three the dual polynomials (in the eigenvalue S{n)) are introduced 
by a different solution method (the three term recurrence (13. 2p ) of the same eigenvalue 
problem (12. ID . The duality at many different levels, e.g. x ^ n, rj^x) ^ ^(n), is displayed 
(I3.14p - (l3.18p . Section four explores the underlying symmetry properties of the Hamiltonian, 
shape invariance §4.11 closure relation §4.21 The general formulas of the coefficients An and 
Cn of the three term recurrence are derived. The forms of the sinusoidal coordinate and 
the diagonal elements of the Hamiltonian are determined in §4.31 from the closure relation. 
The dual closure relation is discussed in §4.41 Section five provides the explicit forms of 
various quantities derived in previous sections for the concrete examples of the orthogonal 
polynomials in the orders given in the review [13]. Some orthogonal polynomials not listed 
in the review [I3| are also discussed in some detail in §5.31 Appendix A gives the possible 
forms of the Jacobi matrices (the Hamiltonians) flA.15p - flA.16p as the solutions of the closure 
relation. The results show the boundedness and unboundedness of the Hamiltonians and 
eigenvalues clearly. Appendix B provides the collection of the definitions of basic symbols 
and functions for self-containedness. Throughout this paper we use the parameter q in the 
range < g < 1. 
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2 Hamiltonian 

The starting point is the eigenvalue problem for a hermitian matrix 7i 

Hi) = 8ij. (2.1) 
The rows and columns of Ti are indexed by non-negative integers x and y, either finite 

case (I) x,y = 0,1,2,. ..,N, (2.2) 

or infinite: 

case (II) x,y = 0,1,2,.... (2.3) 

Let us call the hermitian matrix 7i Hamiltonian, since the eigenvalue problem (12.11) can be 
considered as the Schrodinger equation which is a difference equation, instead of differential 
in ordinary quantum mechanics. The Hamiltonian we consider in this paper has a general 
form 

7^ = -v/^O^e^/DO^- v/I^e-^ (2.4) 
in which the two functions B{x) and D{x) are real and positive but vanish at the boundary: 

B{x) > 0, D{x) > 0, D{0) = ; B{N) = for case (I). (2.5) 

In ( 12.41) (9 == ^ is the differential operator and its exponentiation gives finite shifts: 

e^^^{x) = tp{x±l). 

Thus, in fact, the eigenvalue problem (12. ip can be written as a difference equation on non- 
negative integer points: 

{B{x) + D{x))ij{x) - ^/B{x)D{x + l)i;{x + 1) - ^/ B{x - l)D{x) i){x - 1) = Sipix), 

x = 0,l,...,{N),.... (2.6) 

The boundary condition D{0) = is necessary for the term V'( — 1) does not appear, and 
B{N) = is necessary for the term iIj{N + 1) does not appear in the finite dimensional 
matrix case. It is also easy to see that ipi^) does not vanish for any eigenvector 

^(0)^0, (^(A^) ^ for case (I)). (2.7) 
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Let us remark that the (g-)Askey- scheme of hypergeometric orthogonal polynomials with 
absolutely continuous weight functions, for example, the Meixner-PoUaczek, the Wilson and 
the Askey- Wilson polynomials, have already been well understood as the eigenvalue problem 
of hermitian (self-adjoint) operators Ti p!3| [TO] . The generic forms of the Hamiltonian Ti 
are: 

n ^V{x) e-'^^V{x)* + ^V{x)* e'^^V{x) - V{x) - V{x)\ (2.8) 
n ^/V{x)q-'^^/V{xY + ^/V{xYq'^^/V{x)-V{x) - V{xy. (2.9) 

These also give difference Schrodinger equations, but the shift is in the imaginary direction, 
instead of real in the case (12.61) . 

Although the Hamiltonian Ti (12. 4p is presented in a difference operator form, it is in fact 
a real symmetric tridiagonal (Jacobi) matrix: 

H = (Hx^y), 'Hx,y = Ti-y^x, (2.10) 

Hx^y = -^/B{x)Dix + l)Sx+i,y - ^B{x-l)D{x)5x-i,y + {B{x) + D{x))5x^y. (2.11) 

It is well-known that the spectrum of a Jacobi matrix is simple; that is, there is no degeneracy 
of the eigenvalues 

£(0) <8(l) <8{2) < ■ ■ ■ . 

The first step for solving the eigenvalue equation (12.11) is to rewrite the Hamiltonian (12.41) 
in a factorised form: 

n = A^A, (2.12) 
= y/B{x) - ^D{x)e-^, (2.13) 
A= ^B{x) ~e^^D{x), (2.14) 

with the forward (backward) shift operator A {A^) being hermitian conjugate of each other. 
Or in the matrix form, the tridiagonal (Jacobi) matrix (12.111) is decomposed into a product 
of A^ and A 

{A^)x,y = VB{^Sx,y - ^/D{^6x^l,y, (2.15) 

Ax,y = ^^B{x)5x,y - ^D{x + l)5x+i,y, (2.16) 
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with A"^ (A) having non- vanishing entries on the diagonal and one below (above) the diagonal. 
Throughout this paper we adopt the (standard) convention that the (0, 0) element of the 
matrix is at the upper left corner. The factorisation also means that the Hamiltonian (12.41) is 
semi-positive definite. In fact, the ground state (the lowest eigenvector) <f)o{x) is annihilated 
by A and thus it is a zero-mode of the Hamiltonian: 

AMx) = ^ HM^) = 0, ^(0) = 0. (2.17) 

As we will see shortly, 0o(a^)^ is the orthogonality measure for the eigenpoljTiomials (I2.33p . 
This is the same situation as in the ordinary quantum mechanics and also in the above 
mentioned theory (I2.8p - (l2.9p for the Wilson and Askey- Wilson polynomials, etc. [T3| [TO] . 
The above equation (I2.17p . being a two term recurrence relation. 



Mx + 1) / 5(x) 



0o(x) \iD{x + iy 

can be solved elementarily with the boundary (initial) condition 

0o(O) = 1, (2.19) 



M^) = Sli^^y ^ = 1,2,.... (2.20) 

With the standard convention nfc=n * = the expression (12.200 is valid for x = also. For 
the infinite matrix case (II), the requirement of the finite i"^ norm of the eigenvectors 

00 00 X—l „/ N 

x=0 x=0 y=0 ' 

imposes constraints on the asymptotic behaviours of B{x) and D{x). 

Next let us determine the other eigenvectors (the excited states) in a factored form 

^P(x) = (f)o{x)v{x). (2.22) 

The eigenvalue problem (12. ip for ip is rewritten to that for v 

nv{x)=£v{x), (2.23) 

in which Ti is the similarity transformed Hamiltonian in terms of the ground state wavefunc- 
tion (po(yX): 

'H = <Po{x)-^ono<p^{x) 
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= 5(x)(l-e^) + D(x)(l-e~^). (2.24) 
It gives rise to a simple difference equation for v{x) (12.231) . Obviously a constant is a solution 

7^1 = 0. (2.25) 
In the matrix form Ti. is another tridiagonal matrix 

'H = i'Hx,y), 'Hx,y = B{x){5.j:^y - + D{x){5.^^y - 5.^^i^y). (2.26) 

For all the cases discussed in this paper, the similarity transformed Hamiltonian Ti is 
lower triangular with respect to the special basis 

1, r/(x), vix?, ...,r/(x)", (2.27) 

spanned by the sinusoidal coordinate ri{x) fl3\ [TO]: 

7^?7(x)" = £{n)ri{x)^ + lower orders in ri{x). (2.28) 

Here £{n) is the corresponding eigenvalue. This property is preserved under the affine 
transformation of ^^(x), ri{x) aT]{x) + b {a,b: constants). As will be shown shortly, the 
functional form of the sinusoidal coordinate is characterised by the closure relation (I4.3ip . 
fl4.32p up to a multiplicative and an additive constant. We choose the additive constant to 
achieve 

r/(0) = 0. (2.29) 

The lower triangularity (12.281) implies that the eigenvectors can be obtained as polynomials 
in T]{x) up to normalisation: 

HPMx)) = S{n)Pn{r]{x)), 

n = 0,l,2,...,(iV),.... (2.30) 

H(j)n{x) = S{n)(f)n{x), (f)n{x) = 00 (x)P„ (//(o;) ) , 
We choose the simple normalisation of Pn{ri{x)) 

P„(0) = 1, n = 0,l,2,...,(iV) (2.31) 

which is consistent with (12.71) and (I2.29p . Note that this also means 

Po = 1. (2.32) 



Thus the eigenvectors are completely specified and the eigenvalue problem (12.11) is completely 
solved. The orthogonality of the eigenvectors of the Jacobi matrix Ti (12. 4p . (12.111) implies 

J2 0o(x)'Pn(r?(x))P^(r/(x)) = ^ 5n,m, (2.33) 

in which the normalisation constants {dn > 0} are to be calculated. In terms of the nor- 
malised eigenvectors 

^o{x) = do(j)o{x), 0„(x) = dn(l)oix)Pnir]{x)) = (j)o{x)^Pnir]{x)), (2.34) 

"0 

the orthonormality relation reads 

dn(t)o{x)Pn{ri{x)) ■ dm(t)Q{x)Prn{ri{x)) = Sn,m- (2.35) 

X 

This in turn implies the completeness relation 

dn(t)Q{x)Pn{ri{x)) ■ dn(l)o{y)Pn{ri{y)) = Sx,y, (2.36) 

n 

which states the simplest duality that the row eigenvectors are also orthogonal. When written 
slightly differently 

J2dlPnivix))Pniv{y)) = lAz^2^^^y^ (2-37) 

„ H^O{x) 

it displays the duality with (I2.33P explicitly. 

We have established that Pn{v) is an orthogonal polynomial of a discrete variable with 
the discrete measure {4>o{x)'^} given explicitly as (I2.20p . (I2.2ip . They satisfy the three term 
recurrence relation 

r]Pnir]) = AnPn+M + BnPniv) + CnPn-M (2-38) 

with coefficients An, Bn and C„ (Co = 0). The boundary condition P„(0) = 1 (I2.3ip implies 

Bn = -{An + Cn). (2.39) 

In section m we will derive the general formulas of the coefficients An and C„ (I4.49p . (I4.50p 
from the input functions B{x) and D{x) with the help of various quantities implied by the 
closure relation. This will provide complete specification of the polynomial Pn{v)- 

In the next section, we will show that the explicit form of the normalisation constants 
{dn} (I2.33P can be easily obtained as the components of the ground state wavefunction of 
the dual polynomial as seen clearly from (I2.33P and (I2.37p . 
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3 Dual Polynomials 



The dual polynomial is an important concept in the theory of orthogonal polynomials of a 
discrete variable |2l [3l |4l [Sj |14j. Here we show that the dual polynomial arises naturally as 
the solution of the original eigenvalue problem (]2.ip or (12.23^ obtained in a different way. 
Let us rewrite the similarity transformed eigenvalue problem T-Cv{x) = S v{x) ( ]2.23p into an 
explicit matrix form with the change of the notation v{x) ^{Qo, Qi, • • • , Qx, ■ ■ •) 

J2'^-^yQy = ^Q-^ x = 0,l,...,{N),.... (3.1) 

y 

Because of the tridiagonality of Ti, it is in fact a three term recurrence relation for Qx as a 
polynomial in S: 

8Qx{S) = B{x){Qx{S) - Qx+i{S)) + D{x){Qx{S) - Q.-i(^)), 

x = 0,l,...,(iV),.... (3.2) 

Starting with the boundary (initial) condition 

Qo = 1, (3.3) 

which is consistent with (12.71) . we determine Qx{£) as a degree x polynomial in £. It is easy 
to see 

g,(0) = l, x = 0,l,...,(iV),.... (3.4) 

When £ is replaced by the actual value of the n-th. eigenvalue £{n) (I2.28P in Qx{£), we 
obtain the explicit form of the eigenvector 

Y,'>ixMS{n)) = £{n)Qx{£{n)), x = 0, 1, . . . , (iV), .... (3.5) 
y 

In the finite dimensional case (I), {Qq{£), . . . ,Qn{£)} are determined by (13. 2p for x = 
0, . . . , N — 1. The last equation 

£Qr,{£) = D{N){Qn{£) - Qn-i{S)), (3.6) 

is the degree + 1 algebraic equation (characteristic equation) for the determination of 
all the eigenvalues {£{n)}. This is a familiar situation encountered in quasi-exactly solv- 
able quantum mechanics [151 IIHl [13 IHj- the next section we will provide further two 
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independent algebraic methods for the determination of the eigenvalues {S{n)} based on 
the shape-invariance p3] and the exact Heisenberg operator solution [10] , both of which are 
applicable to the finite (I) as well as the infinite dimensional case (II). 

We now have two expressions (polynomials) for the eigenvectors of the problem (12.11) 
belonging to the eigenvalue S{n); Pn{j]{x)) and Qx{£{n)). Due to the simplicity of the 
spectrum of the Jacobi matrix, they must be equal up to a multiplicative factor a„, 

Pn{r]{x)) = anQA£{n)), x = 0, 1, . . . , (iV), . . . , (3.7) 

which turns out to be unity because of the boundary (initial) condition at x = (I2.29p . 
dMH), (O (or at n = ^M). dM])); 

Pnivm = Pn{0) = 1 = Qo{S{n)), n = 0, 1, . . . , (AT), . . . , (3.8) 

Poivix)) = 1 = QAO) = QxiSm, x = 0, 1, . . . , (iV), .... (3.9) 

We have established that two polynomials, {P„(?7)} and its dual polynomial {Qxi^)}, 
coincides at the integer lattice points: 

Pn{vix))=Q.{S{n)), n = 0,l,...,(iV),..., x = 0, 1, . . . , (iV), . . . . (3.10) 

The completeness relation (12.371) 

J2dlPrMx))PrMy)) = Y.dlQ.{£{n))Qy{£{n)) = -1- (3.11) 

n n VOy ) 

is now understood as the orthogonality relation of the dual polynomial Qx{£), and the 
previous normalisation constant d\ is now the orthogonality measure. 

The real symmetric (hermitian) matrix 'Hx,y (12.111) can be expressed in terms of the 
complete set of the eigenvalues and the corresponding normalised eigenvectors 

'Hx,y = ^S{n)(l)n{x)(t)n{y), (3.12) 
n 

(j)n{x) = dn<pQ{x)Pn{r]{x)) = dn(t>o{x)Q ^{S (u)) . (3.13) 

The very fact that it is tridiagonal can be easily verified by using the difference equation for 
the polynomial Pn{r]{x)) or the three term recurrence relation for Qx{^{n)). 
Here is the list of the dual correspondence: 

X <-> n, (3-14) 
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B{x) ^ -An, 
D{x) ^ -Cn, 



r^(0) = 0^£(0) =0, 



(3.15) 
(3.16) 
(3.17) 



(3.18) 



00 (0) rfo' 



In the last expression, we inserted 0o(O) = 1 fl2.19p for symmetry. It should be remarked that 
B[x) and D{x) govern the difference equation for the polynomial Pn{fl)i the solution of which 
requires the knowledge of the sinusoidal coordinate ri{x). The same quantities B{x) and D{x) 
specify the three term recurrence of the dual polynomial Qx{£) without the knowledge of 
the spectrum £{n). It is required for them to be the eigenvectors of the eigenvalue problem 
fl2.ip . Likewise, An and C„ in (12.381) specify the polynomial Pnijl) without the knowledge 
of the sinusoidal coordinate. As for the dual polynomial Qx{£{n)), An and C„ provide the 
difference equation (in ra), the solution of which needs the explicit form of £{n). Let us stress 
that it is the eigenvalue problem (12.11) with the specific Hamiltonian (12. 4p that determines 
the polynomial Pnijl) ^-^d its dual Qx{£), the spectrum £{n), the sinusoidal coordinate rj{x) 
and the orthogonality measures (t)o{xy and d\. 

4 Shape Invariance and Closure Relation 

Like the Hamiltonians for the (g-)Askey- scheme of hypergeometric orthogonal polynomials 
with continuous measure, all the Hamiltonians for the orthogonal polynomials of a discrete 
variable discussed in this paper are endowed with two types of symmetries; shape-invariance 
and closure relation [10]. The former leads to the determination of the entire spec- 
trum (14. 5 p and the corresponding eigenvectors (14. 6p . In other words, it guarantees the exact 
solvability in the Schrodinger picture. The closure relation (I4.3ip . on the other hand, is 
essential for the exact Heisenberg operator solution (14.340 . which in turn gives the annihi- 
lation/creation operators (I4.37p . They help to determine the coefficients An and C„ (I4.49p . 
(I4.50p in the three term recurrence relation of the polynomial Pn{v) &s well as the algebraic 
method (14.420 to determine the eigenvalues {£{n)}. In §4.31 we show that the requirement 
of the closure relation determines the sinusoidal coordinates 'r]{x) and B{x) and D{x). In 
Appendix A, the possible forms of the Hamiltonians are determined (lA.15p - (]A.16p . 
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4.1 Shape Invariance 

Shape invariance dictates the parameter dependence of the Hamiltonian. Let us denote the 
set of parameters symbohcally by A: 

n{x-X) =A{x-\yA{x;X), or 7i:(A) = ^(A)U(A). (4.1) 

Shape invariance simply means 

Ai\)A{Xy = kA{X + S)^A{X + S)+ S{1 ; A), (4.2) 

in which k is a positive constant and £^(1;A) is the eigenvalue of the first excited state 
£^(1) > with the explicit parameter dependence. That is, the original Hamiltonian 7Y(A) 
and the associated Hamiltonian ^(A)^(A)^ in Crum's [7] sense (or in the factorisation 
method [8], or in the so-called super-symmetric quantum mechanics [9]) have the same 
shape with a multiplicative factor The necessary and sufficient condition for (14.21) is 

^/B{x + 1 ; X)D{x + 1 ; A) = K^/B{x ; X + 6)D{x + l;X + 6), (4.3) 
B{x ■,X) + D{x + l;X) = k{B{x ; X + S) + D{x ; X + S)) + £{1 ; A). (4.4) 

For the finite dimensional case (I), the new (associated) Hamiltonian 7i(A + 5) = ^(A -|- 
6yA{X + 6) is of dimensional, and fOj) holds for < x < - 1. 

It is rather straightforward to extract the necessary information, as done explicitly for 
the (g-)Askey-scheme of hypergeometric polynomials with continuous measure [13] . The 
spectrum is simply generated by £^(1, A): 

n-l 

^(n;A) = ^/t^^(l; A + S(5), (4.5) 

and the corresponding eigenvectors are generated from the known form of the ground state 
eigenvector 0o (I2.2UI) together with the multiple action of the successive A^ operator: 

0„(x ; A) oc ^(A)U(A + (5)U(A + 25^ ■■■AiX+{n- l)S)Uo{x ; A + n(5), (4.6) 

which is related to a Rodrigues type formula. It should be stressed that shape invariance, 
(14. 5 p guarantees the monotonous increasing property (within the valid parameter range) 

^ The continuous g-Hermite polynomial [31 [H] has Hamiltonian (|2.9p with V{x) = 1/ ((1 — z^)(l — qz'^)), 
z — e'^, which has no shiftable parameter. The above shape invariance relation gives a simple realisation of 
the g-oscillator algebra AJ^ = q~^J<^ A + £\, Z\ — q^^ — 1, [19j . 
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of S{n). For all the polynomials discussed in this paper, the spectrum is either linear or 
quadratic in n or g-quadratic in n: 



(i): 


£{n) 


= n, 




(ii): 


S{n) 


= m{n + a), 


e 


iii): 


S{n) 


= 1 - q\ 




iv): 


S{n) 


= q'" - 1, 




(v): 


S{n) 


= 6(g-"-l)(l-ag"), 


e 



1 for a > —1, 
-1 for a < -N, 



1 for a < g ^, 
— 1 for a > . 



(4.7) 

(4.8) 

(4.9) 
(4.10) 

(4.11) 



As we will show shortly in the next subsection, the requirement of closure relation restricts 
the structure of the sinusoidal coordinate ri{x), which has the same form as above with n 
replaced by x, (14. 72 p - (14. 76 p reflecting the duality x ^ n, rj{x) £{n). 

Forward and Backward Shift Operators The counterparts of the A and A'^ operators 
in the polynomial space, called the forward and backward shift operators, play an important 
role in the theory of orthogonal polynomials [3]. Since the shift in x is closely related 
with the parameter shift, we introduce an auxiliary function ^{x) = ip{x ; A) > 0, v^(0) = 1 
to absorb the effects of the parameter shift in the ground state wave function: 



^{x-X) 



def /S(0;A) 0o(a;;A + (5) 



B(x\X) 



(x ; A) 



(4.12) 



In the finite dimensional case (I), B{N ; A) = but it is canceled by the zero in (poiN ; A + (5) 
and ip is well defined. With the help of fl2.18p . it can be rewritten as 



V5(x; A) 



B{0;\) 0o(x;A + ^) 
D{x + l;\) 0o(a: + 1; A)' 



(4.13) 



Explicitly it reads 



ipix;X) 



\B{x-X)^J-^ B{y-X) D{y + l;X + Sy 



(4.14) 



For the Hamiltonians (12.81) . (12.90 with a continuous variable x, the corresponding v?(x) is 
given in [10]. 
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Let us introduce the forward shift operator T = ^(A), and the backward shift operator 
B = B{X) as 

n{\) = S(A)J^(A), (4.15) 
^(A) = ^/B{0;X) 0o(x ; A + 6)-' o ^(A) o ; A), (4.16) 
B{X) ^-1-— 0o(x ; A)-i o ^(A)t o 0o(x ; A + 5). (4.17) 

With the help of fl4.12p - fl4.13p we obtain expressions 

^(A) = 5(0 ; \)^{x ■ X)-\l - e^), (4.18) 
B{X) = [B{x ■,X)~D{x; A)e-^)^(x ; A). (4.19) 

With 5(0; A) appearing in many places, these formulas look rather contrived. However, 
their action on the polynomial Pn{v) = Pnivi^ ! A) ; A) is rather simple: 

^(A)P„(r/(x ; A) ; A) = /„(A)P„_i(r7(x ; X + 6) ; X + d), (4.20) 
B{X)PMx ; A + 5) ; A + (5) = 6„(A)P„+i(r/(a; ; A) ; A), (4.21) 

in which /^(A) and 6n(A) are constants satisfying £{n;X) = 6„_i(A)/„(A). In fact, for all 
the examples given in ^ they take simple values, 

/„(A) =^(n;A), 6„(A) = 1. (4.22) 

Moreover, the auxiliary function <f{x) is related to the sinusoidal coordinate r]{x) as 

. + - (4.23) 



which gives 

r^(l) 



'^"^ (4.24) 



The action of A, A'^ on the eigenfunction 0„(x) is 



^(A)0„(x ; A) = -y^== fn{X)(t>n-i{x ; A + 5), (4.25) 
-d(0 ; A) 



^(A)Vn(a; ; A + (5) = /B(o7A) 6„(A)0„+i(x ; A). (4.26) 

In terms of the forward and backward shift operators, the shape invariance condition 
721) reads 

J'^(A)i3(A) = kB{X + 5)T{X + 5)+ £{l ; A) (4.27) 
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The necessary and sufficient condition for (14.271) are 

B{x + 1 ; X)<^{x + 1 ; A) = kB{x ; A + d)'^{x ; A), (4.28) 
D{x ; X)ip{x - 1 ; A) = kD{x ; A + d)ip{x ; A), (4.29) 
B{x ; A) + L>(x + 1 ; A) = k{B{x ; X + S) + D{x ; X + 6)) + £{l ; A), (4.30) 

wliicli look simpler tlian (14. 3 p and (14. 4p thanks to the presence of (p{x). 
4.2 Closure Relation 

The other symmetry of the Hamiltonians for the orthogonal polynomials is called the closure 
relation. It is realised if the 7i = 7i(A) and the sinusoidal coordinate rj = r]{x ; X) satisfy 

[n, [n,v] ] = vRo{n) + [h,v] Ri{n) + R-i{n), (4.31) 

in which Ri{z) = Ri{z]X) is a polynomial in z. The naming is due to the fact that 
ri{x) + R^i{H) / Ro{H) undergoes a sinusoidal motion with frequency ^^/Ro(H) in the classical 
dynamics limit [10]. By similarity transformation in terms of (polx ; A), (I4.3ip reads 

[n, [n, v] ]=vRo{n) + [n, v] Ri{n) + R-i{n). (4.32) 

The l.h.s. consists of the operators e^^, e^, 1, e~^, e~^^ and Ti contains e^^. Thus Ri can be 
parametrised as 

R^{z) = ri^'z + rf\ Ro{z) = rf\^ + r^^^ z + r^°\ R.i{z) = r^^lz^ + r^^.z + /"J. (4.33) 

The coefficients rp^ depend on the overall normalisation of the Hamiltonian and the sinu- 
soidal coordinate: oc ^(O)^--'', rj^'^ oc ^(O)^-^ ri^'J oc 7/(1)5(0)2"^'. 

Heisenberg operator solution and Annihilation and Creation operators The clo- 
sure relation (14.311) enables us to express any multiple commutator 

[H, [H,--- ,[n,v{x)]--]] 

as a linear combination of the operators ri{x) and [Ti, ri{x)] with coefficients depending on the 
Hamiltonian Ti. only. The exact Heisenberg operator solution for the sinusoidal coordinate 
ri{x) [lOj is given by {t is the time variable): 

e'''^r]{x)e-''^ = a(+)e*"+(^)* + a^^^e*"-^^)* - R.i{H)Ro{n)-\ (4.34) 
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a±{z) = \{Ri{z) ± V^i(z)2 + 4i?o(^)), (4.35) 
Ri{z) = aj^{z) + a^{z), Rq{z) = —a^{z)a^{z)^ (4.36) 
a(±) ±(^[^^^(3,)] _ [r]{x) + R_i{n)Ro{n)-^)a^{n)^[a+{n) - a.{n)Y^ (4.37) 

= ±(a+(7^) - a-{n)Y^ ([n, r]{x)] + a±{n) {i]{x) + i?_i(H)i?o(^)~')) • (4.38) 

The positive and negative frequency parts of the Heisenberg operator solution, a^~^ and 
a^^\ are the annihilation and creation operators, which are hermitian conjugate to each 
other Applying fl4.34p to the eigenvector 0„(x) = (f)o{x)Pn{r]{x)) and using the 

three term recurrence relation for the polynomial Pniv) (12.381) . we obtain [10] 



a(+Vn(a;) = A„0„+i(x), (4.39) 
a(-)</)„(x) =C„0„_i(x). (4.40) 

The constant part in (14.341) is related to 

R^,{S{n))Ro{S{n))-' = -5„ = + C^. (4.41) 

The exact energy eigenvalues {£{n)} can be obtained as the solution of the non-linear equa- 
tion 

a±{£{n)) = £{n±l) - £{n) (4.42) 

starting with £{Q) = [lOj. The above relation simplifies the form of the annihilation- 
creation operators applied to 0„(x) 



a(^Vn(a;) (4.43) 

in) 



.(. + 1)-W-1) ^^^^^ ^ ^'^^^ - ^ ^^^^^^^ ^(^^f^ 



This is the eigenvector version of the raising/lowering operators or the structure relation for 
orthogonal polynomials [221. The polynomial version is obtained by the similarity transfor- 
mation in terms of 0o(3;), aS^'^ 



def 



60 (x) ^ o a(=^) o 00 (a;), 



aWp„(r/(x)) (4.44) 
- .(. + 1)-W-1) (^^'^^^^^ ^ ^'^^^ -'^^ ^ ^^^^^^^ + .(^^1)-%) )^"^^^-^^- 
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Determination of An and C„ The coefficients An and C„ of the three term recurrence 
relation of the eigenpolynomial Pn{v) can be determined elementarily from the above expres- 
sion for Bn (14.411) and the duality. The starting point is the three term recurrence relations 
of the polynomial Pn{v) its dual Qx{£): 

r/(x)P„(r/(x)) = An{Pn+Mx)) - Pnivix))) + Cn{Pn-Mx)) ' Pnivix))) , (4.45) 

S{n)Q,{S{n)) = B{x){Q,{S{n)) - Q,^,{S{n))) + D{x){Q,{S{n)) - Q,^,{S{n))). (4.46) 

By putting x = in (I4.46P and using Qo = 1, D{0) = 0, we obtain 

Pnivil)) = Qi{S{n)) = (4.47) 

in which the ffist equality is due to the duality (13.101) . Next, by putting x = 1 in (14.451) and 
using (14.471) . we obtain 

7]{l){£{n) - B(0)) = An{£{n + 1) - £{n)) + Cn{£{n - 1) - £{n)). (4.48) 

The two equations (I4.4ip and (14.481) for n > 1 give 

S{n+l)~S{n-l) ' ^ ' ' 

^ _ (^H - + ^)) + ^(^) (^H - ^(Q)) 

- S{n-1)-S{n+1) ' ^^-^"^ 

and for n = we obtain from (14.411) 

Ao = R.i{0)Ro{0)-' = r'^'l/rl^\ Co = 0. (4.51) 

They are slightly simplified (ri > 1) in terms of a±: 

_ R.,{S{n)) + v{l){S{n)-B{0))a^{S{n)) 

a4S{n)){a^{S{n))-a4S{n))) ' ^ ' ^ 

^ ^ R-i{S{n)) + v{l){S{n)-B{0))a4S{n)) 

Note that (14.481) for n = gives another important relation 

AoS{l)+B{0)7]{l) =0. (4.54) 
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This relation ensures that the expression for Cn ( I4.53P vanishes for n = 0. Written differently, 
the above relation fl4.54p means an equality 

m-Z^ (4 55) 

^(1) - ,(1) ^^-''^ 

between the two intensive quantities. That is, they are independent of the overall normali- 
sation of the Hamiltonian Ti of the polynomial system and of Ti'^'^^' of its dual. 

Dual Hamiltonian and Determination of dn/do By the duality 03.100 . the three term 
recurrence relation for Pn{i]) (I4.45P translates to the difference equation in n for the dual 
polynomial Qx{S{n)); 

v{x)Q,{S{n)) = A^{Q,{S{n + 1)) - Q,{S{n))) + C^{Q,{S{n - 1)) - Q,{S{n))). (4.56) 

With abuse of notation, it can be written as an eigenvalue equation for the operator Ti'^"'^' 
with eigenvalue 77(0;): 

^dual def _ _ ^^^^ _ ^^^^^^ 

n^'^^'Q^iSin)) = r]ix)Q,i£in)). (4.58) 
The corresponding Hamiltonian matrix is parametrised by n, m = 0, 1, . . . , (A^), . . ., 
H-^"^' = -x/^e^-V^^- v/^e-^"V^^-(A„ + C„), A„ < 0, C„ < 0, (4.59) 



("W )n,m — ~\/ AnCn-\-l 5n+l,m ~ A^-lCn Sn-l,m ~ {An + Cn) Sn,m- (4.60) 

Following the same path that led to (f)o{x) fl2.20p . the orthogonality measure for the polyno- 
mial Pn{f]{x)), we obtain the formula for the orthogonality measure dn for the dual polyno- 
mial Qx{,£{n)): 

dn 



A 



n (4-61) 



m=0 ^"^+1 



do ^ 

It should be stressed that this same dn is the normalisation constant of the polynomial 
Pn{ri{x)) fl2.33p . With the above formula, one only needs to evaluate do explicitly. 

4.3 Determination of r]{x) and B{x) + D{x) 

In the previous subsection, consequences of the closure relation are explored. Here we will 
consider the closure relation as algebraic constraints for unspecified Hamiltonian {i.e. B{x) 
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and D{x)) and as yet undetermined sinusoidal coordinate 77 (x). Under very mild assumptions 
(on top of the general conditions fl2.5p ). rj{^) = 0, ri{x) > rj{y) for x > y, x,y E Z+, 
we determine the possible forms of t]{x) and B{x) + D{x). The characterisation of all the 
systems satisfying the closure relation [i.e. determination of B{x) and D{x) separately) will 
be relegated to Appendix A. 

The necessary and sufficient condition for the closure relation in the polynomial space 
^M) is 

if]{x + 2) - 2r]{x + 1) + r]{x) = rf^r]{x) + /^J + rj^^ [7]{x + 1) - r]{x)) , (4.62) 
ri{x — 2) — 2ri{x — 1) + ri{x) = r^^ri{x) + r^] + r\^^ (^r]{x — 1) — f]{^)) ; (4.63) 
{ri{x + 1) - ri{x)) {B{x + 1) + D{x + 1) - B{x) - D{x)) 

= -(ri^^r/(x) + r^^l) {B{x + 1) + D{x + 1) + B{x) + D{x)) - {r^a\{x) + r^l\) 

~ {ri{x + 1) - ri{x)) (^rS^^ {B{x + 1) + D{x + 1)) + rf ^) , (4.64) 
(r7(x - 1) - ri{x)) {B{x - 1) + /^(x - 1) - B{x) - D{x)) 

= - {r^Q^ri{x) + ri^i) {B{x - 1) + D{x - 1) + 5(x) + D(x)) - (ri^^riix) + r^) 

- {ri{x - 1) - r/(a;)) (^rj^^ (5(s - 1) + D{x - 1)) + rf^) , (4.65) 
- 2(7y(x + 1) - r]{x))B{x)D{x + 1) - 2(77(2; - 1) - r]{x))D{x)B{x - 1) 
= (^S^^??(a;) + r'^-l) (^{B{x) + /^(x))^ + B{x)D{x + 1) + D{x)B{x - 1)) 
+ (rj^^r7(x) + ri\^) {B{x) + D(x)) + rj°^r7(x) + ri°J 

+ ^ (^(r/(x + 1) - r]{x))B{x)D{x + 1) + (r7(x - 1) - r7(2;))D(x)5(a; - 1)) . (4.66) 

The first two equations fl4.62p and fl4.63p . two above and below the diagonal of fl4.32p . 
determine r]{x). Then the next two equations (14.640 and (14.651) . one above and below the 
diagonal, give B{x) + D{x) as a function of 7]{x). In Appendix A, we show that the last 
equation (14.660 . the diagonal part, determines B{x) and D{x) separately. 
By subtracting (14.620 from (14.630 (with x replaced by a; + 2), we obtain 

= {rf^ - r^o^){ri{x + 2) - ri{x)) , (4.67) 

which gives a general constraint on the parameters r^^^ = By subtracting (14.640 from 
(I4.65P (with X replaced by x + 1) and using the above constraint, we obtain 

= (2ri°^ - rj^^) {ri{x + 1) - ri{x)) . (4.68) 
19 



With the two constraints 



(2) _ ^(1) ^(1) _ 9^(0) 

1 5 '0 ~ ^' 1 5 



Tq = r 



(4.69) 



the first equation fl4.62p is reduced to a three term recurrence relation for ri{x): 

ri{x + 2) - (2 + rP)r/(x + 1) + = ri^J, (4.70) 

which can be solved by the input ri{0) = and ^^(1). The solutions are classified by the roots 
of the 'characteristic equation': 

- (2 + rf^)t + 1 = 0. 

In order to have two positive roots (necessary for the condition ri{x) > ri{y) for x > y), 
r^^ > is needed. For the case r[^^ = 0, we have a quadratic solution: 



rj{x) = ^r'^^lx^x - 1 + (for r'^^i ^ 0); r]{x) = r]{l)x (for = 0). (4.71) 

For the rest r^^ > 0, we have a g-quadratic solution. It should be noted that the overall 
multiplicative factor is irrelevant. After choosing 77(1) and other parameters properly, we 
arrive at 



,(2) 



(2) 



(i) : r]{x 

ii) : r]{x 

rj{x 
rj{x 

(v) : rj{x 



[III 
(iv 



X, 

e'x{x + d), 
q-^ - 1, 

eV'-l)(l-cig^), 



1 for > -1, 
-1 for d < -N, 



(4.72) 

(4.73) 

(4.74) 
(4.75) 

(4.76) 



1 foT d < q ^, 
— 1 for d > q~^ . 

Here as usual g is < g < 1. Similar conclusion was reached in a different context by 
Atakishiyev-Rahman-Suslov [2T]. This result can be considered as a difference equation ver- 
sion of Bochner's theorem [221 Ell • fact, the original setting of Bochner's theorem [22] 
is rather limited and its generalisation within the framework of ordinary quantum mechanics 
is given in [10], which contains various sinusoidal coordinates ri{x) = x, x^, cosx, sinhx, e~^, 
etc. 

In these cases the parameters in fl4.69p read 

(i) : rj^^ = 0, ri^J = 0, (4.77) 
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(ii) : ^ 


= 0, ri1 = 


2e', 






(4.78) 


(1) 

111) : r{ 


= (?~^ - 


)^ 


^(2) _ 
'-1 — 


-{q-\ -qhf, 


(4.79) 


iv) : r^^ 


= {q'^ - 


)^ 


^(2) _ 
'-1 ~ 


(q-k -qhf, 


(4.80) 


(v) : ^ 




)^ 


^(2) _ 
'-1 — 


(g-|_gi)2e'(l + rf), 


(4.81) 



and ri{x ± 1) — ri{x) are given by 



(i): 


rj{x ± 1) 


— Tji^x) 


= ±1, 


(4. 


82) 


(ii): 


rj{x ± 1) 


— rj{x) 


= ±e'{2x±l + d), 


(4. 


83) 


iii) : 


rj{x ± 1) 


— rj{x) 


= (1 - 


(4. 


84) 


iv) : 


rj{x ± 1) 


— rj{x) 


= (g^i-l)g-^ 


(4. 


85) 


(v): 


Tji^X ± 1) 


— ri{x) 


= e'(g^i-l)g-"(l-c/g2-±i). 


(4. 


86) 



Solving fl4.64p and (14.651) is somewhat involved. With the constraints (I4.69P incorporated, 
( KM reads, 

{r]{x + 1) - r]{x)) (a^+i - a^) 

= -(r[^^?7(x) + ri^])(a^+i + a^) - {2rf^r]{x) + ri\^) - {r]{x + 1) - //(a;)) (rf ^a^+i + rf^), 

which can be rewritten as 

(ri{x + 1) — ri{x) + r^^ri{x + 1) + ri^])a^+i + (— ?7(x + 1) + ri{x) + r^^ri{x) + r^^l^ax 
= -{ri{x + l)+ri{x))rf -r^y^. (4.87) 

Here we have introduced the abbreviation 

a^ = B{x)+D{x). (4.88) 

With the help of (1470|) . fOTD can be simplified to 

(r/(x + 2) - 77(3; + l))a^+i + (77(2; - 1) - r/(x))a^ = -(r/(a; + 1) + ?7(a;))rf^ - ri\\ (4.89) 

By multiplying r]{x) — r]{x + 1) to fl4.89p . we obtain a two term recurrence relation for 
= B{x) + D{x): 

{r]{x + 2) - r]{x + 1)) (rj^x) - r]{x + l))ar,+i - {r]{x + 1) - t]{x)) {r]{x - 1) - r]{x))(^x 
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= rf ^ {ri{x + 1)^ - ri{xf) + ri\^ {ri{x + 1) - ri{x)) , (4.90) 
which can be solved with the boundary (initial) condition 

ao = 5(0). (4.91) 



The result is that the diagonal component of H and H, =^ B{x) + D{x), multiplied by 
(?7(x + 1) — ?7(x)) (?7(x — 1) — r7(x)) is a quadratic polynomial in ri{x): 

(ri{x + 1) — ri{x)^ (ri{x — 1) — 'ri{x)^arc = rf'\{x)'^ + r^}lri{x) + ri{l)ri{—l)ao, (4.92) 

in which we have introduced a constant //(—I) 

Vi-l)=r^'l-va), (4.93) 

such that (14.701) is formally satisfied. This is a very important characterisation of B{x) and 
D{x) as we will see in the concrete examples in section [5l It should be remarked that the 
factor (^r]{x + 1) — vi^)) {vi^ ~ ^) ~ vi^)) is also a quadratic polynomial in r]{x), 

— (ri{x + 1) — vi^)) {v{^ ~ 1) ~ V{^)) = r^^ri{x)^ + 2ri^^?7(a;) — ri{l)ri{—l), (4.94) 

for all the cases listed in fl4.72l) - (l4.76l) . as can be easily verified from fl4.82l) - (l4.86l) . Therefore 
we obtain 

rf\{xf + r^y^r]{x) + r]{l)r]{-l)B{0) 
r^^r]{xy + 2r^^lr){x) — ri{l)ri{—l) 
The analysis of the remaining equation fl4.66p . which is the diagonal part of the closure 
relation fl4.32p is relegated to Appendix A. 

4.4 Dual Closure Relation 

Since the sinusoidal coordinate plays the central role in the theory of orthogonal polynomials, 
understanding its general properties is important. Let us now consider the 'dual' closure 
relation by interchanging t] ^ Ti (<^=> £) in the closure relation (14.311) : 

[r^, [r/, n]] = n i?;J"^'(r/) + [r/, H] i?f ^^(r^) + (r^), (4.96) 

in which Rf^'^^{T]) are polynomials in rj. By similarity transformation in terms of (f)o{x), we 
obtain 

[r/, [r^, H]]=n i?o^-'(r/) + [r?, H] i?f ^'(r/) + (r/), (4.97) 



_ '1 'IK-^J "T ' -I'/y.^; -r '/v-l;'/v-x;-^vw; 
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which is again a relationship among tridiagonal and diagonal matrices. It is equivalent to 
the following set of three equations: 

{7]{x + 1) - 7]{x)Y = R^''''\v{x + 1)) - {7]{x + 1) - r7(x))i?f "'(r/(a; + 1)), (4.98) 
{r]{x - 1) - 7]{x)Y = R^''''\r]{x - 1)) - - 1) - 7]{x))Rf''''\r]{x - 1)), (4.99) 
= {B{x) + D{x))R^''^\r]{x)) + i?Y(^(a;)). (4.100) 



These imply 



Rf'''\r]{x)) = {r]{x + 1) - r]{x)) + {r]{x - 1) - r]{x)) , (4.101) 
R^'"'\ri{x)) = - {ri{x + 1) - ri{x)) (r/(x - 1) - ri{x)), (4.102) 
i?Y(r^(x)) = -{B{x) + D{x))R^''^\r]{x)). (4.103) 



^dual 




(1) , (2) 


^dual 




= r?z^ + 2ri'; 


^dual 







These expressions are the dual formulas of £{n)^ (14.361) . (I4.4ip and (14.421) . The right hand 
sides of ( I4.101l) - (l4.103p are in fact polynomials in rj{x): 

(4.104) 

ri{l)r]{-l), (4.105) 
7(l)r7(-l)fi(0). (4.106) 

The first two relations are consequences of (14.701) (see also (I4.94p ). and the last relation is 
nothing but the main result (14.921) of subsection 14.31 The coefficients of Rf^^^{z) defined 
like as (14.331) satisfy the same constraints as (I4.69p . The number of independent parameters 
in Rf^^^{z) is six, just the same as in Ri{z) with (I4.69p . Among them, as mentioned after 
(I4.33p . the overall normalisation of Ti and ri{x) are immaterial. This reduces the number 
of essential parameters to four, which corresponds to the most generic cases of the Racah 
§5.1.11 and the g- Racah polynomials §5.1.51 



5 Concrete Examples 

Here we will discuss the properties of various orthogonal polynomials in our scheme. We stick 
to the standard notation as far as possible; the inevitable deviation is caused by symmetry 
consideration and our universal normalisation of polynomials (12.3 ip . (13. 4p . When we use 
non-standard parametrisation (Racah, g-Racah, Hahn, g-Hahn and their duals) we adopt 
different symbols so that no confusion will ensue. We are well aware of the shortcomings 
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of deviating from the standard notation; for example, the connection to other existing con- 
cepts/quantities are blurred. However, we strongly believe that presenting the entire theory 
of orthogonal polynomials of a discrete variable in the hermitian matrix (or operator) form 
has enough merits of elucidating their common structure to compensate the shortcomings. 
The orthogonal polynomials whose orthogonality relation is written in terms of a g-integral 
(Jackson integral), e.g. the Big g-Jacobi polynomial, etc. will need different treatment and 
will be discussed elsewhere. We start with the finite dimensional case (I) and then move to 
the infinite dimensional case (II). In each case the Askey-scheme of hypergeometric orthogo- 
nal polynomials will be followed by the g-scheme polynomials, in the order presented in the 
review of Koekoek and Swarttouw [14J, with the number e.g. [KS3.2] attached indicating 
the subsection there. The subsection is named after the polynomial to be discussed in the 
subsection. We explain various quantities and concepts in some detail for the first examples, 
the (g-) Racah polynomials. For the rest, mostly formulas only are given except for those 
needing some attention and/or remarks. As mentioned earlier, a good part of the explicit 
formulas for specific polynomials is already known in various contexts. 

The only input is the Hamiltonian 7i or the two positive functions B{x) and D{x). In 
order to define them, the parameters must be specified. So we start from the parameters and 
their shift properties under shape invariance. The other quantities are defined and derived 
as explained in previous sections. Most basic definitions and notation are recapitulated in 
Appendix B for self-containedness. 

5.1 Finite Dimensional Case (I) 

In this subsection x takes a finite range of values 

x = 0, 1, 2, . . . , iV. 
5.1.1 Racah [KS1.2] (self-dual with different parameters) 

As well-known, the Racah polynomial is the most general hypergeometric orthogonal poly- 
nomial of a discrete variable. All the other (non-g) polynomials are obtained by restriction or 
limiting procedure. The set of parameters A and their shifting unit S and the multiplication 
factor K, (14.21) are: 

X={a,b,c,d), ^ = (1,1,1,1), K = l, (5.1) 
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which are different from the standard ones 7,5). Let us define d, e and e' 



d'=a + b + c — d — 1, 



dcf 



1 for d> -1, 



^ -1 for d< -N, ^ 



dcf 



1 ioT d > -1, 
-1 for d < -N, 



(5.2) 
(5.3) 



and B{x) and D{x) as 

B{x]X) = -e 
D{x;X) = -e 



{x + a){x + b){x + c){x + d) 
{2x + d){2x + I + d) ' 
{x + d — a){x + d — h){x + d — c)x 



{2x - I + d)i2x + d) 
The finiteness condition B{N) = can be achieved by 

a = —N or b = —N or c = —N. 



(5.4) 
(5.5) 



(5.6) 



Since B[x) and D{x) are symmetric in a, b and c, we assume without loss of generality 



-iV, a>b. 



(5.7) 



Then we restrict our argument to the following parameter ranges in which B{x) and D{x) 
are positive: 



(i>0, a>A^ + rf, 0<6<l + (i 

d < -2N, a>0, N + d<b<l~N 

d>0, 0<a<l + d, b<l-N 

d < -2N, N + d<a<l-N, b<l + d -- 

The energy eigenvalue and the sinusoidal coordinate are 



(e,e' 



(1,1), 
(1,-1), 



-1,1), 



-1,-1). 



£{n ; A) = en{n + d), rj{x ; A) = ex{x + d). 



(5.8) 
(5.9) 
(5.10) 
(5.11) 



(5.12) 



The polynomial is 



P„(r/(x ; A) ; A) = 4F3 



-n, n + d, —X, X + d 
a, 6, c 



Rn{e'rj{x ;A);a — l,rf — a, c — l,(i — c) 



(5.13) 



in which 4F3 is the standard hypergeometric series flB.3P and i?„ is the standard notation 
for the Racah polynomial [H]. The orthogonality measure (or the squared ground state 
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wavefunction) 4>o{xy and the normalisation constants (12.331) are obtained from (12.201) 
and (I4.6ip (with An and C„ given below), 

'^^^'^'^^ - {l + d-aA + d-hA+d-cA)^^r' ^^-^^^ 



dn{\) 



2 {a,b,c,d)n 2n + d 



[1 + d — a,l + d — b,l + d — c, 1)„ d 



-1)^(1 +d-a,l + d-b,l + d-c) 



(5.15) 



{d + l)Nid+l) 2N 

in which (a)„ is the Pochhammer symbol fIB.ip . The format for c?„(A)^, throughout this 
section, consists of two parts separated by a x symbol: d^ = (d^/dl) x d^ reflecting the 
duality (I3.18p . The second part d^ satisfies the relation J2x^o(-^y ~ ^/^o- The coefficients 
of the three term recurrence (12.381) An and Cn for the polynomial Pn{v) obtained from the 
closure relation (I4.49p - (l4.50p are: 

^ ^ ^, {n + a){n + b){n + c){n + d) 
{2n + d){2n + 1 + d) ' 
^ ^ ^, (n + J- a)(n + (i - &)(n + J - c)n ^ 
(2n-l + d)(2n + (i) 

The three functions -Ri, Rq and -R_i appearing in the closure relation (14.310 are: 

Ri{z;\) = 2e, (5.18) 
Ro{z ; X) = Aez + d^ - 1, (5.19) 
R.i{z ; A) = 2eV + ee' {2{ab + be + ca) - {1 + d){l + d))z + e'abc{d- 1). (5.20) 

Under the parameter shift A ^ A + 5 the functions B{x) and D{x) behave: 

B{x;X + S) = B{x + l;X)^±^, D{x ; X + 5) = D{x ; X)^^^. (5.21) 

From this and (I4.14p . (p{x) is 

^ ^ 2x + d + 1 , ^ 

^ix;X)= . (5.22) 

The constants /n(A) and bn{X) appearing in (I4.20l) - (l4.2ip are 

fn{X)=S{n;X), bn{X) = l. (5.23) 

By comparing B{x), D{x) with An and C„, (poi^)'^ with d^, it is clear that the theory is 
self-dual with the parameter correspondence (a, b, c, d, e, e') (a, b, c, e', e). 
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5.1.2 Hahn [KS1.5] (dual to dual Hahn ^5.1.31) 

The parameters (a, 6) are very slightly different from the standard ones [a, (3) (a + 1 = 
a, f3 + 1 = b) for the Hahn polynomial. For obvious reasons, we adopt the same parameters 
for the Hahn and dual Hahn polynomials: 



A=(a,6,iV), 5 = (1,1,-1), K = l, 

_ r 1 for a, 6 > 0, 
^ ~ I -1 for a,6 < 1 - A^, 

B{x ; A) = e{x + a){N - x), D{x ; A) = ex{h + N -x 

£{n ; A) = en{n + a + b — 1), rj{x ; A) = x, 

-n, n + a + 6 — 1, —x 



P4v{x;X);X) = 3F2[ ^ _^ ' | l) = Qn(r/(x ; A) ; a - 1, 6 - 1, iV), 

(j)o{x;Xy 



^2 N\ (a), (6)jv-x 



x\{N-x)\ (6)jv 

dnixy 



2 A^! {a)n{2n + a + b-l){a + b)N {b)N 



n\{N-n)\ (6)„(n + a + 6- l)jv+i {a + b)N 

{n + a){n + a + b - 1){N - n) 
{2n-l+a + b){2n + a + b) ' 
n{n + b - l){n + a + b + N - 1) 



AniX) : 

'^"^'^^^ {2n-2 + a + b){2n-l + a + b)' 
Ri{z;X) = 2e, 

Rq[z ; A) = 4ez + (a + 6 - 2)(a + 6), 
R-i{z ■ A) = -t{2N - a + b)z - a{a + b~ 2)N, 
B{x;X + S) = B{x + 1; A), D{x;X + d) = D{x ; A), 
ip{x-X) = l, fn{X) = £{n-X), bn{X) = l. 



5.24) 
5.25) 

5.26) 
5.27) 
5.28) 

5.29) 

5.30) 

5.31) 

5.32) 

5.33) 
5.34) 
5.35) 
5.36) 
5.37) 



Obviously the Hahn and dual Hahn polynomials are dual to each other. The standard 
parameters for the latter are (7,5). 

5.1.3 dual Hahn [KS1.6] (dual to ^5.1.21) 

The parameters (a, b) are very slightly different from the standard ones (7, 5) (7 + 1 = 
a, (5 + 1 = 6) for the dual Hahn polynomial: 

A=(a,6,A^), (5 = (1,0,-1), /t = l, (5.38) 

^={ \ f M (5-39) 

I -1 for a, 6 < 1 - iV, ^ ' 
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{x + a){x + a + b - 1){N - x) 
(2x - 1 + a + b){2x + a + b) ' 
x{x + b - l){x + a + b + N - 1) 



B{x ; A) 

^(a; ; ^) - ^2a; - 2 + a + 6)(2x - 1 + a + 6) ' 
£{n ; A) = n, ?7(x ; A) = ex{x + a + 6 — 1), 

-n, X + a + 6 — 1, —X 
a, -A^ 



P„(r7(x; A); A) = 3^2 



(5.40) 

(5.41) 
(5.42) 

1 ) = Rn{eT]{x ; A) ; a - 1, 6 - 1, iV), (5.43) 



0o(a;; A) 
rfn(A)2 = 



{a)^{2x + a + b- l){a + b) 



N 



x\{N-x)\ {b)^{x + a + b-l)N+i 

N\ (a)„ {b)N-n (&)iV 

— X 



n\{N-n)\ {b)N (a + 6)^' 

An{X) = -e{n + a){N - n), C„(A) = -en{b + N -n), 

Ri{z;X)=0, Ro{z;X) = l, R^i{z ; X) = 2ez^ - e{2N - a + b)z - eaN, 



(5.44) 

(5.45) 
(5.46) 
(5.47) 



B{x;X + S) = B{x + l; A) ^^,1 ^ | t ^ D{x-X + d) = D{x ; A) ^""^ ^ ^ t ^ (5-4^ 



2x + a + 6 

y?(x; A) = , /„(A) = ^(n; A), 6„(A) = 1. 

a + 



2x + a + 6 



(5.49) 



With the present parametrisation, the following duality fl5.28p ^ (15.431) . (15.271) ^ (I5.42p . 
(15:26|) ^ dMSD, fl53T|) - (l5:32|) ^ (ICTjl - dCTj) . dOOD ^ ([525]), flCT]) ^ (ICTD . is obvious. 

5.1.4 Krawtchouk [KSl.lO] (self-dual) 



A=(p,iV), 5 = (0,-1), K = l, 
< p < 1, 

5(x; A) =p(A^-x), D{x]X) = {l-p)x, 



£{n ; A) = n, //(x ; A) = x, 

-iV 



P,(r/(x; A); A) = 2^1 
AT! 



-1 



ir„(r/(x;A);p,A^), 



^0 



>; A) 



x! (A^ — x)! V 1 — p 



N\ 



h(t^)"'^(^-^) 



iV 



n! (Ar-n)!Vl-j9. 
A„(A) = -p(Ar - n), C„(A) = -(1 - p)n, 

Pi(^;A) = 0, Po(^;A) = l, R^i{z-X) = {2p-l)z-pN, 

E(x ; A + (5) = 5(x + 1 ; A), L>(x ; A + 5) = D{x ; A), 

<^(x;A) = l, /„(A) =^(n;A), 6„(A) = 1. 



5.50) 
5.51) 
5.52) 
5.53) 
5.54) 

5.55) 

5.56) 
5.57) 
5.58) 
5.59) 
5.60) 



28 



5.1.5 g-Racah [KS3.2] (self-dual with different parameters) 

This is the first example of the g-scheme of the orthogonal polynomials. Among them the q- 
Racah polynomial is the most general. The set of parameters A is different from the standard 
one 7,5) in the same manner as for the Racah polynomial. The shifting unit 6 and 

the multiplication factor k (14.21) are: 

= (a,6,c,rf), (5 = (1,1,1,1), K = q-^. (5.61) 

Here stands for q^^^'^^^---) = (q^^, q^^, . . .). They are shifted multiplicatively. By definition 
the parameter q is not shifted. Let us introduce d, e, e' 



d = abcd-^q-\ (5.62) 

1 for < q~^, / _ f 1 for (i < q~^, 
-1 for J>g-^, ^ ~ I -1 for > g"^. 



(5.63) 



The functions B(x) and D(x) are 



The finiteness condition B{N) = can be realised by one of the following choices: 

a = or b = q^^ or c = q~^ . (5.66) 

Thanks to the symmetry in a, h and c, we choose for simplicity 

c = g-^, < a < 6, < rf, (5.67) 

and consider the following parameter ranges in which B{x) and D{x) are positive: 

Q<d<l, Q<a<q^d, qd<h<l ^ (e, e') = (1, 1), (5.68) 

d > g-2^, < a < 1, q^-^ <b<q^d (e, e') = (1, -1), (5.69) 

< < 1, grf < a < 1, 6 > g^-^ ^ (e, e') = (-1, 1), (5.70) 

d > g-2^, q^-^ <a< q^d, b > qd (e, e') = (-1, -1). (5.71) 

In case some of the parameters a, b and d are non-positive, the situation is slightly compli- 
cated but can be treated in a similar way. The energy eigenvalue, the sinusoidal coordinate 
and the g-polynomial are: 

S{n ; A) = e(g~" - 1)(1 - Jg"), r?(x ; A) = e'{q-' - 1)(1 - rfg^), (5.72) 
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P„(r7(x ; A); A) =403 



x 

^ q;q 



a, b, c 

= Rn{e'{l + d + r]{x ; A)) ; aq^^, da~^, cq~^, dc~~^), (5.73) 

in which 4^3 is the basic hypergeometric series (lB.4p . The orthogonahty measure (poix)"^ and 
the normahsation constants d'^ ( 12.33^ are obtained from ( ]2.20p and ( 14.6ip (with An and C„ 
given below), 

; A)^ = Mc,d;q). ^ 1-d^ ^^^^^^ 

{a~^dq,b~^dq,c~^dq,q;q)xd^ 1 — a 
^ ^^^2 _ {a,b,c,d;q)n 1 - Jg^" 



a ^dq,b ^dq,c ^dq,q;q)nd"- 1 — d 

(-l)^(a-Mg, b-'dq, c'^dq ; q)N rf^g^^(^+^) 



(5.75) 



{dq;q)N{dq;q)2N 

in which (a; g)n is the g-Pochhammer symbol flB.2p . The coefficients of the three term 
recurrence (12.381) An and C„ for the polynomial Pn{v) obtained from the closure relation 
f H:i9|l - fl430|) are: 

^„(;,)^,,(1-«^")(1:''^")(1-<)(1-''V'), (5,75) 

. ,/l~a-^dqnil-b-^mil-c-^^^^^^^ 
[1 - dq^^-^){l - dq^") 

The three functions Ri, Rq and appearing in the closure relation (14.311) are: 

R,{z;X) = {q-^ -q"2)''z', z' = z + e{l + d), (5.78) 
Roiz;X) = (g-^ -gi)2(z'2 - (g"^ +g^)2j), (5.79) 
R_i{z]\) = (g^5 -qhf(^e'{l + d)z'^-ee'{a + b + c + d + d+{ab + bc + ca)q~^)z' 

+ e'((l-a)(l-6)(l-c)(l-rfg-i) 

+ {a + b + c-l-dd+ {ab + bc + ca)q^'^)){l + d)). (5.80) 



For the shifted parameters the functions B and D are: 

B{x;X + d)=B{x + l;X)^—^^^, D{x ; X + 5) = q'D{x ; X) (5.81) 
From this and fl4.14p . (p{x) is 

^(x;A) = ^^^^^. (5.82) 
30 



The constants /n(A) and bn{X) appearing in fl4.20l) - fl4.21l) are 

UX)=S{n;\), 6„(A) = 1. (5.83) 

By comparing B{x), D{x) with An and C„, (t)i){xY with d^, it is evident that the theory is 
self-dual with the parameter correspondence (a, 6, c, d, e, e') (a, 6, c, d, e', e). 

5.1.6 g-Hahn [KS3.6] (dual to gSHHD 

Here and in the next subsection (dual g-Hahn) we will consider, for simplicity, the positive 
parameter range of a and 6, which are very slightly different from the standard ones (a,/5) 
{aq = a, f3q = b) for the Hahn and dual Hahn polynomials. The standard parameters for 
the latter are (7,5). These are the multiplicative counterparts of the parameter change for 
the Hahn polynomials. The parameter range could well be enlarged. 



g^=(a,6,g^), 5 =(1,1,-1), = g-\ 

1 for < a < 1, < 6 < 1, 
-1 for a > q^-^, b > g^"^, 

B{x ; A) = 6(1 - ag^)(g^-^ - 1), D{x ; A) = eaq-\l - g^)(g^-^ - b), 
S{n ; A) = 6(g-" - 1)(1 - abq^-'), ^{x ; A) = q'^ - 1, 

= <5n(l + ^?(a;; A) ;ag"\6g"\A^|g), 

{q;q)N {a;q)x{b;q)N-x 



:x;xy 



{<i;(l)x{q;q)N-x {b;q)Na^ 



^ ^ iq;q)N {a,abq-^;q)n 1 - a&g^""^ ^ {b;q)Na^ 

{q;q)n{q;q)N-n {abq^,b;q)na'' 1 - abq^^ {ab;q)N 
(^qn-N _ i)(i _ ag")(i _ a6g"-i) 



C„(A) = -aq 



(1 - a6g2n-i)(i _ a5g2n) 

(1 - g'^)(l - a6g"+^-i)(l - 6g"-i) 
(1 - a6g2n-2)(i _ abq^"^-^) 



Ri{z ; A) = (g-2 - ga)^/, 2' "i^' z + e{l + abq 
R,{z ■ A) = (g-5 - g^)2(^'2 _ ^^(1 ^ ^-1)2)^ 

i?_i(z ; A) = (g-^ - g^)2(^^'2 _ ^^^^ _^ ^^-i) _^ ^ ag-i)g-^)e/ 

+ a(l + g-i) ((a - l)6g-i + (1 + 6g-i)g-^ 
B{x;X + d) = B{x + 1;\), D{x;X + d) = q^D{x ; A), 
¥?(x ; A) = g""^, /„(A) = ^(n; A), 6„(A) = 1. 



5.84) 

5.85) 

5.86) 
5.87) 

5.88) 
5.89) 
5.90) 
5.91) 

5.92) 

5.93) 
5.94) 

5.95) 
5.96) 
5.97) 
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Obviously the g-Hahn and dual g-Hahn are dual to each other. 



5.1.7 dual g-Hahn [KS3.7] (dual to ^5.1.61) 

For obvious reasons, we adopt the same parameters (a, b) for the g-Hahn and dual g-Hahn 
polynomials: 



a. a / 



{q ; q)x {q ; q)N~x {abq^, h\q)xa'' 1 - ahq- 
^ /^^2 ^ (g;g)jv (a;g)n(fe;g)jv-« ^ {b;q)Na^ 

{q;q)n{q;q)N-n {b;q)Na'^ {ab;q)N 
A4X) = -e(l - ag")(g"-^ - 1), ^(A) = -eaq-\l - g")(g"-^ - b), 

i?i(;2;A) = (g-5-gi)V, z' = z + 1, 

Ro{z-X) = {q-^~q^fz'\ 

R.i{z ; A) = (g-^ - g^)'(e(l + abq-')z''' - e(a(l + bq-') + (1 + ag-i)g-^) 

+ ea(l + g^i)g-^ 



l-a6g2^+2 l-a6g2^-2 



i?(a;;A + (5) = i?(x + l;A)^— D{x ; X + S) = q'D{x ; X) ^ _ ^^^2, 

V{^'^^) = —. 7—^ fnW=^{n;X), bn{X) = l. 

1 — ab 

5.1.8 quantum g-Krawtchouk [KS3.14] 

g^=(p,g^), (5 = (1,-1), «: = g, p > g"^, 

5(x ; A) = p-ig^(g^-^ - 1), D{x ; A) = (1 - g^)(l - p-\^-^-^), 



(5.99) 



g^ = (a,6,g^), 5 = (1,0,-1), k = q-\ (5.9 

1 for < a < 1, < 6 < 1, 
-1 for a > gi-^, b > g^^^, 

, (g-^-l)(l-ag-)(l-a&g-^) 

^"^'^ (l-a6g2-i)(l-aV) 

nr. - - -.,x-iv-i (i - - afeg^-+^-^)(i - bq^-') 

^ ^ (l-a6g2-2)(l-a6g2-i) 

^(n;A) =g-"-l, r]{x;X) = e(g-^' - 1)(1 - afeg^-^) 

'g~", abq^~^, q~^ 
a, q- 

= i?„(e(l + abq-^ + r]{x ; A)) ; aq~\ bq-\ N\q), 
2_ {q\q)N {a,abq-^ ■,q)^ 1 - afog^^-^ 
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£(n;A) = l-g", r/(x ; A) = g-^ - 1, 



6o(x ; A 



{q;q)x{q]q)N~x {p^^q ^■,q)x' 
c^n(A) = ^ , . _„ — ^ X [p q ; q)N 



{q;q)n{q;q)N-n {p ^q ";g)r 

A„(A) = -j9-ig-"-^-i(l - g^-"), C„(A) = -(g"" - 1)(1 -p-^g"'^), 
i?i(z;A) = (g-^-g^)V, 
i?o(^;A) = (g-^-g^)V^ 

; A) = (g-5 - gi)2(^'2 + + p + g-^-i)/ + + g-i)), 

; A + 5) = g-^5(a; + 1 ; A), D(a; ; A + 5) = qD{x ; A), 
<^(x ; A) = g""", /n(A) =^(n; A), 6„(A) = 1. 

The dual quantum g-Krawtchouk polynomial with the same parameter and 

B{x ; A) = p-ig-^-^-i(l - g^"-), D{x ; A) = (g-^ - 1)(1 - p-'q'^, 
S{n;X) = q-^-l, r]{x ; X) = 1 - q^ , 

AnW = -p^^g"(g"-^ - 1), a(A) = -(1 - g")(l -p^^g"-^-!), 



5.114) 
5.115) 

5.116) 

5.117) 
5.118) 
5.119) 
5.120) 
5.121) 
5.122) 
5.123) 



5.124) 
5.125) 
5.126) 



has not been reported in Koekoek-Swarttouw [H]. It is interesting to note that the functions 
B{x) and D{x) in fl5.124p are related to those of the affine g-Krawtchouk polynomial §5.1.101 
with the change of variable x ^ N — x, B ^ D: 

B{x) = D'^'^"" g-Krawtchouk _ ^ ^(^^^ ^ ^affine g-Krawtchouk^^ _ (5.127) 

with the change of the parameter p^^ffino g-Krawtchouk _^ p-'^q-^-^. 

5.1.9 g-Krawtchouk [KS3.15] and dual g-Krawtchouk [KS3.17] 

g^ = (p,g^), (5 = (2,-l), k = q-\ p > 0, (5.128) 

B{x;\) =q'=-^ -1, D(x; A) =p(l -g^), (5.129) 

£(n;A) = (g-"-l)(l+pg"), r/(x ; A) = g-^ - 1, (5.130) 

P„(r^(x ; A) ; A) = 302 \q;q)= + v{x ; X) ■,P, N ; q), (5.131) 
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{Q;Q)x{q;q)N-x 

Ri(z;\) ^(q~^ -q^fz', z'^z + l-p, (5.135) 

i?o(^;A) = (g-i -g^)2(/2+p(g-i +gi)2), (5.136) 

; A) = {q--2 - q^f{z" + {p - q-'')z' + p{l + q-'){l - g"^)), (5.137) 

S(x;A + 6) = S(x + l;A), D{x ; X + d) ^ q'^D{x ; X), (5.138) 

(/?(x;A) = g-^ /„(A) = £:(n;A), 6„(A) = 1. (5.139) 

The dual g-Krawtchouk polynomial in this parametrisation has 

^^^■M r(x-A)-z..^-^-^ (i-g-)(i+pg-+^) . . 

^^'^'^^-(l+pg2.)(i+p^2.+i)' ^(^>A)-pg (l+pg2.-i)(i+^^2.)' (5.140) 

£:(n;A) = g-"-l, 77(2; ; A) = (g"^ - 1)(1 +pg^), (5.141) 

^n(A) = -(?"-^-l), C„(A) = -p(l-g"). (5.142) 

These are to be compared with the standard parametrisation of the dual g-Krawtchouk 
[KS3.17], in which the parameter c is identified as 

c=-pq^. (5.143) 

dual g-Krawtchouk [KS3.17] in the standard parametrisation 

(c,g^), (5= (0,-1), K^q-\ c<0, (5.144) 

(5.145) 

(5.146) 

£(n;A)=g-"-l, 77(0; ; A) = (g^^ - 1)(1 - cg^"^), (5.147) 

— ti ■ — X X — TV 

P,(77(a;;A);A)=302(^ ' q ; q) = K^{1 + cq-'' + rj{x ; X) ;c, N\q), (5.148) 

v g , u / 

iq;q)N (cg-^ ; g). g^--|-(-+i) 1 - cg^-^^ 

' " (g;g).(g;g)iv-. (cg;g).(-c)^ ' ^'"'''^ 

rfn(A)^ = , [';'\ (-c)-g^"("-^) X ^^i^, (5.150) 

(g ; g)n(g ; g)jv-n (c ^;g)jv 
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^ (l-cg2^-^)(l-cg2^+i-^)' 

Dfr • A) = rn^^-2iv-^ (1 - g")(l - eg") 
^ ^ (l-cg2-i-^)(l-cg2-^)' 



^n(A) = -(g"-^ - 1), C„(A) = cg-^(l - g"), (5.151) 

R^{z;X) = {q-^ - q^fz', z' = z + 1, (5.152) 

i?o(-z;A) = (g-5 _gi)V^ (5.153) 

R.i{z ; A) = (g-^ - gi)2((l + cg-^)^^ - (1 + c)g-V), (5.154) 

5(a:;A + 5) = i?(a; + l;A)^— ^^5^:^, D{x ; X + S) = q'D{x ; X) (5.155) 

—X 1— A'^ a:; 

(^(x;A)= ^ lIcV/ ' /n(A)=£(n;A), MA) = 1- (5.156) 



5.1.10 affine g-Krawtchouk [KS3.16] (self-dual) 



g^=(p,g^), 5 = (1,-1), K^q-\ < p < q-\ 

B{x ; A) = (g^-^ - 1)(1 - pq"^^), D{x ; A) = pg^-^(l - g^), 

£:(n;A) = g-"-l, 77(0: ; A) = g"^ - 1, 

P„(77(a:;A);A) = 302(^ ° g ; g) = i^f (1 + r;(a: ; A) ;p, TV ; g), 

\ pg, g / 

(g ; q)N {pq ; g)x 



^0 



:x;xy 



{Q;(l)x{q;q)N-x ipqY 



^"^^^ "(g;g)„(^;g).-. (m)" ^'''^ ' 

A„(A) = -(g-^ - 1)(1 -pg«+i), Cn{X) = -pg"-^(l - g"), 
i?i(^;A) = (g-^-g^)V, z'^^'z + l, 
i?o(^;A) = (g-^-g^)V^ 

i?_i(^ ; A) = (g-^ - g^)2 {z" - {pq + (1 + p)g^^)/ + p(l + g)g-^) , 
B{x;X + S) = B{x + l;X), D{x;X + S) = q^D{x ; A), 
ip(x ; X) ^ q-^ , fn(X) ^£(n;X), 6„(A) = 1. 

5.2 Infinite Dimensional Case (II) 

In this subsection x takes an infinite range of values 

x = 0,l,2,...,. 

In contrast to the finite dimensional case, the structure of the polynomials is severely con- 
strained by the asymptotic forms of the functions B{x) and D{x), which are determined in 
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Appendix A. It is easy to see for ri{x) = x (i) (14.721) . a quadratic spectrum £{n) ~ n{n + a) is 
not possible. In order this to happen in fl2.28p . B{x) and D{x) must have the same coefficient 
(3 with the different sign for the leading quadratic term, B{x) ~ D{x) ~ — which 
violates the positivity of B{x) and D{x). It is also easy to see rj{x) = e'x{x + d) (ii) (14.731) 
is not possible. If the leading power of B{x) and D{x) is quadratic, the coefficient must be 
the same with the same sign, B{x) ~ D{x) ~ /3 > 0. Then for large n (I2.28P gives 
a negative leading term 

Hr]{xY ~ -l3n^7]{xY + ■ ■ ■ , (5.169) 

which cannot give a normalisable eigenvector for a positive definite Hamiltonian. If the 
leading power of B{x) and D{x) is linear, the results in Appendix A (1A.12P and (lA.14p tell 
that r^i^ = and B{x) and D{x) have the opposite sign leading terms 

Bix) = —— — , , , ~ ; -V + lower order, (5.170) 

^ ' ri{l)B{0){ri{x+l)~ri{x-l)) ' ^ ' 

Dix) = — — , , , ~ ^ :v + lower order. (5.171) 

^ ^ t]{1)B{0){t]{x + 1) ~ r]{x - 1)) ^ ' 

This simply contradicts the positivity of B{x) and D{x). If the leading power of B{x) and 
D{x) is a constant or a negative power in x, then eigenpolynomials in ri{x) of 7i do not exist. 
Thus we have two self-dual polynomials in x, the Meixner §5.2.11 and Charlier §5.2.2[ 

As for the g-polynomials, the consequences of the asymptotic behaviours of B{x) and 
D{x) are easy to see. If these functions are bounded as in the cases of the sinusoidal coordi- 
nates r7(x) = q~^ — l (iv) (14. 75 p . and ri{x) = e' {q^^ — — dq^) (v) (I4.76p . the Hamiltonians 
are bounded and eigenvalues cannot take unbounded forms S{n) = g^" — 1, (g^" — 1)(1 — ag") 
and vice versa. As mentioned in Appendix A, there is no self-dual g-polynomial of infinite 
dimension. 

5.2.1 Meixner [KS1.9] (self-dual) 

This is the ffist example of the self-dual polynomial (I5.175p . which is symmetric under 
interchange x n, with S{n) = n and ri{x) = x (I5.174p . Compare (I5.173p ^ (l5.177p . 
flSTTel) : 

A = (/5, c), (5 = (1, 0), K = l, p>0, < c < 1, (5.172) 
B{x;X) = j^{x + p), D{x;X) = j^x, (5.173) 
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£{n ; A) = n, rj{x ; A) = x, 

-n, —X 



P„(r/(x;A);A) =2^1 



l_c-i) =M„(r7(x;A);/?,c), 
X (1 - cY 



AJX) 



1-c 



1-c 



n, 



1 + c 



i?i(z;A) = 0, i?o(^;A) = l, i?_i(2;A) = -- ^ 

1 — c 

B{x]\ + 5) = B{x + l]X), D{x;X + 6)= D{x ; A), 
<^(x;A) = l, /„(A) =^(n;A), 6„(A) = 1. 



5.2.2 Charlier [KS1.12] (self-dual) 

A = a, (5 = 0, K, = 1, a > 
B{x ; A) = a, D(x ; A) = x, 
£{n ] X) = n, rj{x ; A) = x, 
P„(r7(x;A);A) = 2Fo^"''' ""^ 

^0 



-1 



Cn{r]{x ; A) ; a) 



x; A)2 = ^, rf„(A)2 = ^ X e 
x! 



An{X) = -a, C„(A) = -n, 

Ri{z ; A) = 0, Ro{z ; A) = 1, R-i{z ; X) = —z — a, 
B{x;X + 6) = B{x+l;X), D{x;X + d) = D{x ; A), 
¥.(x;A) = l, /„(A) =^(n;A), 6„(A) = 1. 



5.2.3 little g-Jacobi [KS3.12] 



The universal normalisation (12.311) differs from the standard normalisation of the little q- 
Jacobi polynomial pn as shown explicitly in (15.1931) . This does not affect the normalisation 
measure 0o(3;)^ (I5.194p but the normalisation constants d\ (I5.195P and the coefficients of the 
three term recurrence An and C„ (I5.196I) - (I5.197I) are different from the standard ones: 



q^=[a,b), 5 = (1,1), K = g \ 0<a<g \ 6<g \ 

S(x;A) = a(g-^-6g), D(x ; A) = g"^' - 1, 

S{n ; A) = (g-" - 1)(1 - a6g"+i), r^(x ; A) = 1 - q\ 



(5.190) 
(5.191) 
(5.192) 
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P„(?7(x; A) ; A) = (-a) g 2 v ^ 20i 

(6g ; g)„ V aq 



q ; 



:-a)-V^"("+^^|^Pn(l - r/(a; ; A) ; a, 



U'n[^) — 7 ^ 7 ^ 



(g, ag ; g)„ 1 - abq {ahq^ ; q)^ 



2n+i {l-bq--^'){l-abq--^': 



(1 -a6g2n+l)(l_a5g2n+2) ' 

^ (l-g")(l-ag") 

^ (1 -a6g2n)(l_a5g2n+l) ' 

; A) = (g^^ _ gs)^^;', 2;' =^ 2: + 1 + a6g, 
i?o(^;A) = (g-^-g5)2(/2_^5(l + g)2)^ 

i?_i(2;A) = {q~^ -q^f{-z'^ + {l + a)z' -a{l + q){l-bq)), 
B{x;X + d) = q^B{x + 1 ; A), ; A + 5) = D(x ; A), 

ip{x;X) = q'', fniX) = £in;X), 6„(A) = 1. 

The dual little g-Jacobi polynomial with the same parameter and 



(1 -a6g2^+i)(l-a6g2^+2)' ' ' ' (1 - a6g2^)(l - afeg^^+i^ 
; A) = 1 - g", ; A) = (g-^ - 1)(1 - a6g^+^), 

Ar^iX) = -a(g-" - 6g), C„(A) = -(g-" - 1), 



5.193) 
5.194) 
5.195) 
5.196) 

5.197) 

5.198) 
5.199) 
5.200) 
5.201) 
5.202) 

5.203) 
5.204) 
5.205) 



has not been reported in Koekoek-Swarttouw [T3]. The dual little g-Jacobi polynomial 
was introduced by Atakishiyev and Klimyk [2^, which has the same three term recurrence 
relation and the orthogonality measure as above. The difference equation was not mentioned 
there. 

5.2.4 g-Meixner [KS3.13] 

g^ = (6,c), 5 =(1,-1), K = q, 0<b<q-\ c> 0, (5.206) 
5(x;A) =cg^(l-6g^+^), /^(a; ; A) = (1 - g^)(l + 6cg^), (5.207) 
£(n;A) = l-g", r7(x ; A) = g"^ - 1, (5.208) 

Pn{7]{x ; A) ; A) = 201 ^ ' I g ; -c-^g"+^) = M„(l + r]{x ; X) ■,b,c; q), (5.209) 
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M^^-AY- , <(Af= , .t^^, (5.210) 

{q,-bcq;q)^ [q, -c-^q ; q)n (-c;g)oo 

A„(A) = -eg— - ^(A) = -(g-" - 1)(1 + eg""), (5.211) 

R,{z;X) = {q'^ -q^fz', z' = z - 1, (5.212) 

i^o(^;A) = (g-^-g5)V^ (5.213) 

i?„i(2;A) = (g"5 _g5)2(2'2 + (i_c_5c)z'-c(l + g-^)), (5.214) 

B{x;X + d) = q-^B{x + l;X), D{x ; X + 6) = D{x ; X), (5.215) 

</.(x ; A) = g-", /„(A) =^(n;A), 6„(A) = 1. (5.216) 

The dual g-Meixner polynomial with the same parameter and 

5(x;A) = cg""-^(g-"-6g), D(x ; A) = (g~" - 1)(1 + eg""), (5.217) 
^(n; A) = g-" - 1, r7(x ; A) = 1 - g^ (5.218) 
A„(A) = -cg"(l - 6g"+i), C„(A) = -(1 - g")(l + 6cg"), (5.219) 

has not been reported in Koekoek-Swarttouw [H]. This can be considered another g- version 
of the Meixner polynomial §5.2.11 with the sinusoidal coordinate r]{x) = 1 — q^. 

5.2.5 little g-Laguerre/Wall [KS3.20] (dual to ^5.2.61) 

The universal normalisation (12.311) differs from the standard normalisation of the little g- 
Laguerre/Wall polynomial p„ as shown explicitly in (15.2231) . As other examples, df^, An and 
Cn are different from the standard ones: 



q^ = a, S = 1, K = q ^, < a < q ^ , 
i?(x ; A) = ag~^', L)(x ; A) = g"^' - 1, 
£{n;X) = q-^-l, ^^(x ; A) = 1 - g^ 

Pn(^?(a;;A); A) = 20o(^ g ; a'^g'^J = (a"^g~" ; g)„p„(l - 77(0; ; A) ; a|g), 

1 r ^^2 ((^QT , /^^2 / n 

(/)o(x;A) =- — , dn{X) =- — x(ag;g)oo, 

{q;q)x (g,ag;g)„ 

A4X) = -aq^''^\ C„(A) = -(1 - g")(l - ag"), 
Ri{z-X) = {q-\-q\fz\ z' = z + l, 
R,{z-X) = {q-'^-q^fz'\ 

R-i{z;X) = (g-^ -g5)2(-;z'2 + (l + a)/-a(l + g)), 
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5.220) 
5.221) 
5.222) 

5.223) 

5.224) 

5.225) 
5.226) 
5.227) 
5.228) 



B{x;X + S) = q^B{x + l;X), D{x ; \ + 6) = D{x ; X), (5.229) 
y.(x;A) = g^ fniX)=£in;X), 6„(A) = 1. (5.230) 

The dual little g-Laguerre/Wall polynomial with the same parameter and 

B{x;X) =aq^''+\ D(,t; A) = (1 - g^)(l - ag^), (5.231) 
£(n;A) = l-g", r]{x;X) = q-'' - 1, (5.232) 
A„(A) = -ag-", a(A) = -(g"" - 1), (5.233) 

is Al-Salam Carlitz II polynomial as seen below. This was reported by Atakishiyev and 
Klimyk [25]. 

5.2.6 Al-Salam-Carlitz II [KS3.25] (dual to §[EMD 

The universal normalisation (12.311) differs from the standard normalisation of the Al-Salam- 
Carlitz II polynomial Vn""^ as shown explicitly in (15.2371) . As other examples, d^, An and Cn 
are different from the standard ones: 



q^ = a, 6 = 0, K, = q, < a < g \ 

B{x ; A) = ag2^+\ D{x ; A) = (1 - g^)(l - ag^), 

f(n;A) = l-g^ r]{x ■ X) = q-^ - 1, 

P„(r^(x;A);A) = 20o(^ ^ ^ | = (-a"'?^^""'^)" Ki'^^l + ^(^^ ; A) ; g), 

0o(a;;A)^ = - — , c?„(A)^ = ^x(ag;g)oo, 

{q,aq;q)^ {q;q)n 

An{X) = -aq-\ Cn{X) = -(g"" - 1), 

i?i(z;A) = (g-i-g^)V, ^"=^^^-1, 

i^o(^;A) = (g-5-g^)V^ 

R^,{z-X) = {q-"^-q^\z'^ + {l + a)z'), 

B{x;X + 6) = q-^B{x + 1 ; A), L>(a; ; A + 5) = D(s ; A), 

(^(x; A) = g"'^, /„(A) =f(n; A), 6„(A) = 1. 



5.234) 
5.235) 
5.236) 

5.237) 

5.238) 

5.239) 
5.240) 
5.241) 
5.242) 
5.243) 
5.244) 



5.2.7 alternative g-Charlier [KS3.22] 

The universal normalisation (I2.3ip differs from the standard normalisation of the alternative 
g-Charlier polynomial Kn as shown explicitly in (I5.248p . As other examples, An and Cn 
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are different from the standard ones: 



= a, 6 = 2, K = q ^, a > 0, 
B{x;X) = a, ; A) = - 1, 

£{71 ; A) = (g-" - 1)(1 + ag"), r/(x ; A) = 1 - g^ 



P„(r/(x;A);A) = g"%0i( 



q q X 



q;-a-\-''+^ 





= (-ag")-"ir„(l -r/(x; A) ;a;g), 
(x;A)^= "f ^ , ^„(A)^ = " f ^ / "'^7 ^ + X 



' {q;q)n (-ag";g)oo l + a {-aq;q), 

A^[X) = -aq--- \^ , C„(A) 



2nA 



(1 + ag2«)(l + 052^^+1) ' ' ' (l + ag2"-i)(l + ag 

A) = - gi)22;', z' = Z + I - 

Ro{z;X) = {q-l - q-2f{z'^ + a{q-'^ + q^f) , 
R^iiz ; A) = (q-^ - g5)2(_/2 ^ ^/ _ ^ ^ 

5(a; ; A + 5) = ^^^(a; + 1 ; A), D(a; ; A + 5) = D{x ; A), 
V?(x;A)=g'^, /„(A) = £(n ; A), 6„(A) = 1. 

The dual alternative g-Charlier polynomial with the same parameter and 

5(x;A) = ag3^+^-— ^i±^^^— — , D(x;A) ^ " 



(l + ag2^)(l + ag2^+i)' ^ ' ' (1 + ag2^-i)(l + ag^^) ' 
£(n ; A) = 1 - g", r7(x ; A) = (q-^ - 1)(1 + ag^), 

A„(A) = -a, C„(A) = -(g-"-l), 

has not been reported in Koekoek-Swarttouw [14j. The dual alternative g-Charlier poly- 
nomial was introduced by Atakishiyev and Klimyk [25], which has the same three term 
recurrence relation and the orthogonality measure as above. The difference equation was 
not mentioned there. 

5.2.8 g-Charlier [KS3.23] 

g^ = a, S = -1, K, = q, a > 0, (5.259) 
E(x;A) = ag'^, D{x ; X) = I - q"" , (5.260) 
£(n;A) = l-g", r7(x;A) = g-^-l, (5.261) 
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P„(r/(x ; A) ; A) = 201 Y h ' -«"'?"^') = Cn{l + vi^ ; X) ] a; q), (5.262) 

0o(x;A)^= ^ , d^{Xf = - ^ -X- (5.263) 

A„(A) = -aq-'^~\ CniX) = -(g-" - 1)(1 + aq-^), (5.264) 

R^{z;X) = {q-^ -q^fz', z' = z - 1, (5.265) 

i?o(^;A) = (g-^-g^)V2, (5.266) 

i?_i(2;A) = (g-3 -g^)2(y2^(i_a)2'_a(l + g-i))^ (5.267) 

B{x;X + d) = q^^B{x + l;X), D{x ; X + 6) = D{x ; X), (5.268) 

(^(x;A) = g-^ /„(A)=^(n;A), 6„(A) = 1. (5.269) 

The dual g-Charlier polynomial with the same parameter and 

B{x;X) = aq'^''~\ D(x ; A) = (g^^ - 1)(1 + ag"^), (5.270) 

£:(n; A) = g-" - 1, r/(x;A) = l-g^ (5.271) 

A„(A) = -ag", C„(A) = -(1 - g"), (5.272) 



has not been reported in Koekoek-Swarttouw [14]. This can be considered another g- version 
of the Charlier polynomial §5.2.21 with sinusoidal coordinate ri{x) = 1 — q^. It should be 
emphasised that these two different g-versions of the Charlier polynomial have very different 
characters. The g-Charlier has a bounded spectrum and its dual has an unbounded spectrum. 
This is always the case for the infinite dimensional g-polynomials and their duals. 

5.3 Other Polynomials 

Here we will discuss orthogonal polynomials which are not listed in Koekoek-Swarttouw's 
review [1^. The naming of the polynomials is very tentative. 



5.3.1 alternative g-Hahn 

The g-Hahn polynomial §5.1.6l is a g- version of the Hahn polynomial §5.1.2l with the sinusoidal 
coordinate r]{x) = q~^ — 1. The other g- version with r]{x) = 1 — q^ and the same spectrum 
S{n) = e(g^" — 1)(1 — abq"^^) as the g-Hahn polynomial is as follows (the parameter range 
could well be enlarged): 

(a,6,g^), (5 = (1,1,-1), k = q-\ (5.273) 
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1 for < a < 1, < 6 < 1, 
-1 for a > gi-^, b > q^~^ , 

B{x ■ A) = eaq-\l - - b), D{x ; A) = 6(1 - ag^-^)(g-^ - 1), 

S{n;\) = e(g-"-l)(l-a6g"-^), r/(x;A) = 1 - 

P.(.(x;A);A) 



6, q-^ 

( n ■ n\ »r 

l)o{x ; A) 



2 {<1;(1)n a''{a;q)N-x{b;q)cc 



TV 



2 {q;q)N (6,a6g i;g)„a"g"(" 1) 1 - a6g2n 1 (a ; g);v 



{q;q)n{q;q)N-n {a,abq^]q)n l-abq-^ {ab]q) 

^n{^) — —ciq 



N 



(1 - a6g2n-i)(i _ a5g2n) 
(1 - - ag"-^)(l - a6g"+^-^) 

; A) = (g-^ - g5) V, z' = z + e{l + abq~'), 
R,{z ■ A) = (g-^ - g5)2(/2 _ ^5(1 + g-i)2)^ 

R-i{z-\) = {q--- - q^f[-z'^ + [l + aq-^ + a{l + bq-^)q'')ez' 

- a(l + g-i) (1 - 6 + 6(1 + ag-i)g^) 
S(x ; A + 5) = q^B{x + 1 ; A), L)(x ; A + 5) = D(x ; A), 
¥?(x;A) = g"^, /„(A) = ^(n ; A), 6„(A) = 1. 



5.274 

5.275 
5.276 

5.277' 
5.278 
5.279 
5.280 

5.281 

5.282 
5.283 

5.284 
5.285 
5.286 



This is the most generic form of the functions B{x) and D{x) for the sinusoidal coordinate 
ri{x) = 1 — g^ (quadratic polynomials in g~^) as mentioned (1A.8P in Appendix A. It is 
interesting to note that the functions B{x) and D{x) are related to those of the g-Hahn 
polynomial with the change of variable x ^ N — x, B ^ D: 

B{x) = D''-^''^'^{N - x), D{x) = fi'?-H^h'^(A^ - x). (5.287) 

5.3.2 alternative g-Krawtchouk 

The g-Krawtchouk polynomial §5.1.91 is a g-version of the Krawtchouk polynomial §5.1.41 
with the sinusoidal coordinate ri{x) = q~^ — 1. The other g-version with ri{x) = 1 — g^' and 
the same spectrum S{n) = (g"" — 1)(1 + pq") is as follows: 

g^=(p,g^), 5 = (2,-1), K = q-\ p > 0, (5.288) 
B(x;A) =p(l -g^-^), D(x ; A) = g"^ - 1, (5.289) 
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£{n-\) = (g-" - 1)(1 + M"), r/(x ; A) = 1 - 



n „—x 

— 1 x—N 

q;-p 



P„(r/(x; A) ; A) = 301 

(g;g)^(g;g)7v-x 



N 



(1 +J9g2")(l +pg2n+l) ' '^V / (1 + pg2n-l) (1 + p^2n) 

A) = (g^^ - g5)2/^ z"^^z+l-p, 
Rq{z;X) = {q-"2 - qif[z'^ + p{q~^ + q-^ f), 

R.,iz ; A) = (g-^ - q'^ n-z" + (1 - pg^)/ - p(l + g)(l - g^)) , 
B{x;X + d) = q^B{x + 1 ; A), D{x]\ + 6) = D{x]X), 
(^(x;A) = g^\ /„(A) =^(n; A), 6„(A) = 1. 



5.290) 
5.291) 

5.292) 
5.293) 

5.294) 

5.295) 
5.296) 
5.297) 
5.298) 
5.299) 



The functions B{x) and D{x) are related to those of the g-Krawtchouk polynomial with the 
change of variable x ^ N ~ x, B ^ D: 

B{x) = £)9-Krawtchouk^^ _ ^ j^^^^ ^ ^g-Krawtchouk _ _ (5.300) 

5.3.3 alternative afRne g-Krawtchouk (self-dual) 

The affine g-Krawtchouk §5.1.101 is a self-dual polynomial with ri{x) = q~^ — 1 and S{n) = 
g"" — 1. There is another self-dual polynomial with ri{x) = 1 — q^ and S{n) = 1 — g": 

g^=(p,g^), (5 = (1,-1), K = q, p>q-'', (5.301) 
B{x ; A) = (1 - g^-^)(l - p'^-'), D{x ; A) = p-^g^-^(g-^ - 1), (5.302) 
£(n;A) = l-g", r/(x ; A) = 1 - g^', (5.303) 



-71 „ — X 



P„(r^(x;A);A) = 202(^ '\ (5.304) 

V pg, g / 

0o(x;A)^ ^.^.(..1-^)^ (5.305) 
(g ; g)^ (g ; g)7v-x (p^^g ^ ■,q)N-x 

dniXr = , /^rC'f'' pV^"^^-^^ X (pg)-^, (5.306) 
{q;q)n{q;q)N-n [p ^q ^ ■,q)N-n 

A^{X) = -(1 -p-lg-«-l)(l - g^-"), C„(A) = -p-lg^-"(g-" - 1), (5.307) 

R^{z;X) = {q'^ -q^fz', z' = z - 1, (5.308) 
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Ro{z-X) = {q'-^ -q'^fz'^ (5.309) 

R.,{z ; A) = - q-2f (-/^ _ ^^-i^-i ^ ^ p-^)q^)z' - p-\l + g-^)g^) , (5.310) 

5(x;A + 5) =5(x + l;A), D(x ; A + (5) = g"^D(x ; A), (5.311) 

¥.(x;A) = g^ /„(A) = ^(n ; A), 6„(A) = 1. (5.312) 

The functions B{x) and D{x) are related to those of the quantum g-Krawtchouk polynomial 
with the change of variable x ^ N — B ^ D: 

B{x) = Dq"a>^t""^ g-Krawtchouk(^^ _ jj^^^-^ ^ ^quantum g-Krawtchouk _ (5.313) 

6 Summary and Comments 

A unified theory of orthogonal polynomials of a discrete variable is presented through the 
eigenvalue problem of hermitian matrices of finite or infinite dimensions. The hermitian 
matrices (7i) are real and tridiagonal (Jacobi) matrices (12.41) that can be factorised Ti. = 
A'' A and thus are positive semi-definite. The orthogonality measure ^0(2;)^ is obtained as 
a solution of the factored equation A4>o{x) = 0, as in the ordinary quantum mechanics. 
Then the eigenvalue problem can be solved in two different ways; the first as a difference 
equation and the second through the three term recurrence relations. The former gives 
the eigenvector polynomial Pn{i]{x)) and the latter provides the dual polynomial Qx{^{n)) 
satisfying the relation Pn{r]{x)) = Qx{S{n)) through the universal normalisation condition 
-Pn(O) = 1 = (5x(0). Here S{n) is the eigenvalue of H and ri{x) is the sinusoidal coordinate 
satisfying the closure relation [H, [H,ri{x)]] = ri{x)Ro(H) + [H,ri{x)]Ri(H)+R-i{l-C). Thanks 
to the closure relation and the shape invariance, the entire spectrum, the coefficients 
and Cn of the three term recurrence relation of Pn{v)j its normalisation constants d"^ are 
determined algebraically. As byproducts, some as yet unexplored (not reported in the review 
of Koekoek-Swarttouw [13]) dual polynomials with explicit forms of the difference equation, 
three term recurrence, the normalisation measure/constants are presented in §5.3[ On top 
of them we mentioned the dual quantum g-Krawtchouk fl5.124l) - fl5.126p . the dual g-Meixner 
iKm\\ - (Km\f and the dual g-Charlier (Km\f - iK7m . 

After completing the main part of the paper, we became aware of the work of Terwilliger 
on the correspondence between the Leonard pair [27] and a class of orthogonal polyno- 
mials. Although the general setting of the problem and the methods are markedly different. 
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Terwilliger's results have some overlap with the finite dimensional case (I) of the present 
paper. 
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Appendix A: Possible Forms of the Hamiltonians 

Here we continue from §4.31 the algebraic analysis of the closure relation to determine the 
possible function forms of B{x) and D{x). The remaining equation to be solved is fl4.66p . 
which is the diagonal part of the closure relation fl4.32p . 

With (14.691) and (14. 70 p . the equation to be analysed (I4.66P reads 

(r/(x) - r]{x + 2))B{x)D{x + 1) + (r/(x) - - 2))B{x - l)D{x) 

= {r]{x + 1) - 2ri{x) + r/(x - l))al + {2rf^r]{x) + r'Ll)a^ + rj°^r/(x) + ri°J. (A.l) 

Before proceeding, let us stress again that the input of the eigenvalue problem (12. ip is B{x) 
and D{x). The quantities r-"''' {i = 1,0, —1; j = 0, 1,2) are determined by the consistency 
of the closure relation (I4.3ip . In Appendix A, however, we turn the logic and regard r^''^ as 
input data and determine the possible forms of B{x) and D{x) which satisfy (I4.62p - (l4.66p . 
In this approach the positivity of B{x) and D{x) must be verified at the last step, since it 
is not guaranteed automatically. By multiplying ri{x — 1) — ri{x + 1) to (lA.ip . we obtain 

{ri{x + 2) - r]{x)) (//(x + 1) - r]{x - l))B{x)D{x + 1) - (x a; - 1) 

= [{v{x - 1) - - {v{x + 1) - 

+ (ri{x — 1) — r]{x + 1)) (^2rf^'ri{x) + ri^])a^ + [r]{x — 1) — r]{x + 1)) (r'^^r]{x) + r^^l) 
(rf^rjlx + l)r]{x) + r^}lr]{x + 1) + a'^j (^rf^ri{x + l)ri{x) + r^_^lri{x) + a'^j 

(r]{x + 1) - 

— rf\{x + l)ri{x) — r^^\ (r]{x + 1) + vi^)) — (x x — 1), (A. 2) 
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where Oq = ri{l)ri{—l)aQ and oq = B{0). For the second equahty, the main result of §4.3[ 
fl4.92p is used. The above equation is simphfied as 

(ri{x + 2) - ri{x)) {ri{x + 1) - vi^ - l))B{x)D{x + 1) 

(rj°^?7(x + l)ri{x) + r^_^lri{x + 1) + Oq) (^r^^ri{x + l)ri{x) + r^}lri{x) + Oq) 

[r]{x + 1) - ri{x)Y 

- r7(l)ao(r*:] + ?7(l)ao) - rf\{x + l)ri{x) - r^_^l[r]{x + 1) +r]{x) - //(-I)). (A. 3) 



For each case of ri{x), fl4.72p - fl4.76p . we will solve B{x) from the above equation flA.Sp . 
with the result of a^: = B{x) + D{x), (14. 920 . Let us parametrise B{x) as 

(ri{x + 1) — {v{^ + 1) ~ ~ 1)) 

Note that the denominator is the dual counterpart of a+{S{n)) (a+(£^(n)) — a_(£^(n))) in the 
denominator of An in (14.521) . 

We eliminate D{x + 1) from flA.Sp (multiplied by (rj^x + 2) — rjlx + 1)) (r]{x + 1) — ri{x))'^) 
by using (KW2\i : 

- B{x) ({rf'\{x + if + r5?7(x + 1) + a'^) {r]{x + 2) - r]{x)) + B{x + 1) [7]{x + 1) - r?(a;))j 
= {ri{x + 2) -ri{x + l)) 

X (^rf\{x + l)rj{x) + r*L]?7(x + 1) + dg) (r[°-'r7(x + 1)^7(0;) + rf3]r7(x) + Oq) 
- (?7(x + 1) - rj{x)Y {rj{l)aQ[r^3, + '7(l)ao) + rQ\{x + l)r7(x) 

+ ri°J {r]{x + 1) + 77(3;) - r/(-l)))^ . (A.5) 

This equation determines B{x) uniquely from the boundary (initial) value 5(0) = ri{l){ri{l) 
-r^(-l))5(0). 

We have carried out this last step by assuming that B{x) is a polynomial in x, B{x) = 
YljLoPj^'' foi' the cases (i) and (ii) and that B{x) is a Laurent polynomial in q^, B{x) = 
J2j=Ah (^j^^^ f*^^ cases (iii) - (v). The results are as follows, 
case (i), ri{x) = x: 

A quadratic polynomial (M = 2) solution B{x) is obtained and B{x) and D{x) are: 

B{x) = (quadratic in x), D{x) = x(linear in x), (A. 6) 
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which contains the Hahn polynomial §5.1.21 as the most generic case, 
case (ii), ri{x) = e'x{x + d): 



A quartic polynomial (M = 4) solution B{x) is obtained and B{x) and D{x) are: 
, , (x + (i) (cubic in x) x(cubic in x) 



(A.7) 



{2x + d){2x + I + dy ^' {2x + d){2x - I + dy 



which contains the Racah polynomial §5.1.11 as the most generic case, 
case (iii), rj{x) = 1 — 

A quadratic polynomial (Mi = 0, M2 = 2) solution B{x) in is obtained and B{x) and 
D{x) are: 



The most generic one in this category reported in the Koekoek-Swarttouw's review [14] is 
the little g-Jacobi polynomial §5.2.3[ The dual quantum g-Krawtchouk polynomial fl5.124p - 
(15.1261) also belongs to this case. There is a more general polynomial with one more parameter 
as tentatively called alternative g-Hahn polynomial §5.3. 1[ For the infinite dimensional case 
(II), the functions B{x) and D{x) grow exponentially. Therefore the corresponding spectrum 
is unbounded, i.e. type (iv) (14.101) or (v) (14. lip . In other words, there is no self-dual 
ri{x) = 1 — g^, S{n) = 1 — g" polynomial in the infinite dimensional case (II). 
case (iv), ri{x) = q^^ — 1: 

A quadratic polynomial (Mi = —2, M2 = 0) solution B{x) in q~^ is obtained and B{x) and 
D{x) are: 



which contains the g-Hahn polynomial §5.1.61 as the most generic one. In contrast to the 
previous case the functions B{x) and D{x) are all bounded in the the infinite dimensional 
case (II). Therefore the corresponding spectrum is also bounded, i.e. type (iii) (14.91) only. In 
other words, there is no self-dual ri{x) = q'^ — 1, S{n) = g"" — 1 polynomial in the infinite 
dimensional case (II). 
case (v), T]{x) = e'{q~^ — 1)(1 — dq^): 

A quadratic polynomial (Mi = —2, M2 = 2) solution B{x) in g^ and g~^ is obtained and 
B{x) and D{x) are: 



B{x) = (quadratic in g D{x) = {q ^ — l)(linear in g ^). 



(A.8) 





(A.9) 




(1 — dg^') (cubic in g^) 



D{x) 



(1 — g^') (cubic in g^) 



(A. 10) 



(l-c/g2^)(l-rfg2^'+i)' 



{1 - dq^''){l - dq^'^-^)' 
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which contains the most generic case of the g-Racah polynomial §5.1.51 The functions B{x) 
and D{x) are all bounded in the the infinite dimensional case (II). The corresponding spec- 
trum is also bounded, i.e. type (iii) (14.91) only. The case (v) does not have self-dual polyno- 
mial, either. Thus the most generic functional forms of B{x) and D{x) are determined for 
each case of the five different sinusoidal coordinates (i)-(v). 

Determination of polynomials satisfying certain forms of difference equations for given 
ri{x) (quadratic or g-quadratic) has a long history [231 l2ll [271 126]. Bochner's theorem [22] 
on Sturm-Liouville polynomials is a precursor. The present characterisation in terms of the 
closure relation (I4.32p is consistent with the existing ones. 

It is interesting to note that in all the known cases the following relation is satisfied: 

+£.z7^-(z7^y=o. (A.11) 



5(0)2 5(0)7^(1)5(0)2 \r]{l)B{Oy. 

When this condition is met, the solution of (]A.5|) for all the cases (i)-(v) has a simple 
expression in terms of ri{x): 

B{x) = -rS°^r/(x)r/(x + 1) + ^(^)^(q) ^(^) (^(^ + 1) " vi^)) 

- r^^lriix) + ri{l)B{0){ri{x + I) - ri{x) - ri{-l)) . (A.12) 

For the other function D{x) we parametrise 

D^j,-^ = -^(^^ (A.13) 

(^ri{x — 1) — ~ 1) ^ + 1)) 



The corresponding expression for D[x) is 

M 



D{x) = -rf''^r]{x)r]{x - 1) + -^j^JJ^^^vix) {v{x - 1) - vi^)) 



- r%{x) + r/(l)5(0) {r]{x - 1) - r]{x) - //(-I)) . (A.14) 



The above relation flA.lip is a consequence of fl4.54p . which can be regarded as a con- 
sistency condition. From (gSSD, and fH3T|) we have £{1) = a+(^(0)) = |(rj°^ + 
^(0)2 _|_ ^^(0) -J _ ^(0)^^(0)^ Substituting these into (14.541) (multiplied by rj^^ /£{!)), 
we obtain 



r^l + ^(v/-!°^' + 4rr-rf^)^(l)5(0) = 0, 
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namely 

1 ^ / (0)2 , ^ (0) 1 (0) 

Squaring this gives flA.lip . 



By using (14.541) . the expressions for B{x) and D{x), (1A.4P with (]A.12p and (1A.13P with 



flA.14p . are cast in a similar form as An (14.521) and C„ (I4.53P : 

(r/(x)) +g(l)(r/(a;) + Ao)<-'(r/(x)) 

Bix) = r, ; 7 ; ^-—^ r , A. 15 

af^'(r7(x))(af - al^^i(rj(x))) 

iJ X = r, ; 7 ; ^— ; r . A. 16 

a'^^''\'q{x)){a'^^^\ri{x)) - af^\'q{x))) 

Here R'^^^^{z) and a'^^'^{z) are functions appearing in the dual closure relation §4.4[ The 
functions are given in (I4.1U4|) - (I4.1U6I) and af''\z) are defined by (c.f. f05]) . (02D ) 

^dual(^) dcf 1 j^^dual(^) ^ ^i?f^l(z)2+4/2^nal(^) (A.17) 

which satisfy 

a±^H?7(x)) =r7(x±l) -r7(x). (A.18) 

Similarly An and C„ can be cast into the forms like B{x) in flA.4p with (IA.12P and D{x) in 
with ^KAM . 

Appendix B: Some definitions related to the hypergeo- 
metric and Q'-hypergeometric functions 

For self-containedness we collect several definitions related to the (g-)hypergeometric func- 
tions |[T3] . 

o Pochhammer symbol {a)n '■ 

{a)„ = ff (a + A; - 1) = a(a + 1) ■ ■ ■ (a + n - 1) = (B.l) 
-LJ- r(a) 



In 

k=l 



o g-Pochhammer symbol (a ; g)„ : 

(a ; g)n = - «^'"') = (1 - a)(l - ag) ■ ■ ■ (1 - ag^^). (B.2) 

fc=i 

o hypergeometric series ^Fs '■ 

oo , \ r> 

def (Ctl, ■ ■ ■ ,Clr)nZ' 



ai, ■ ■ ■ , Or 

bi, ■■■ ,bs 



z 



n=0 
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where (oi, ■ ■ ■ , ar)n = IVj=li(^j)n = {ai)n ■ ■ ■ {ar)n- 

o g-hypergeometric series (the basic hypergeometric series) r4>s '■ 

where (ai, ■ ■ ■ , ; g)„ =^ 115=1 ('^i ! = (^i ; ■ ■ ■ (^r ; 
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